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ABSTRACT 



>- ■ We present uniform, detailed spectral analyses of 7-ray burst (GRB) X-ray 



afterglows observed with ASCA, Beppo-SAX, Chandra, and XMM-Newton, and 



CN ' critically evaluate the statistical significances of X-ray emission and absorption 

Q ■ features in these spectra. The sample consists of 21 X-ray afterglow observa- 

^ ■ tions up to and including that of GRB040106 with spectra of sufficient statistical 

quality to allow meaningful line searches, chosen here somewhat arbitrarily to 
P^' be detections with more than 100 total (source plus background) counts. This 

Q ■ sample includes all nine X-ray afterglows with published claims of line detec- 

tions. Moderate resolution spectra are available for 16 of the 21 sources, and 
c^ ■ for the remaining five the Chandra transmission grating spectrometers obtained 

high-resolution data. All of the data are available from the public archive. We 
^ ■ test a simple hypothesis in which the observed spectra are produced by a power- 

j^ ■ law continuum model modified by photoelectric absorption by neutral material 

both in our Galaxy and possibly also local to the burst. As a sample, these 
afterglow spectra are consistent with this relatively simple model. However, 
since the x^ statistic is not sensitive to weak and/or localized fiuctuations, we 
have performed Monte Carlo simulations to search for discrete features and to 
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estimate their significances. Our analysis shows that there are four afterglows 
(GRB011211, GRB030227, GRB021004, and GRB040106) with line-like features 
that are significant at the Sa level. We cautiously note that, in two cases, the 
features are associated with an unusual background feature; in the other two, 
the fractional magnitudes of the lines are small, and comparable to the expected 
level of systematic uncertainty in the spectral response. In addition, none of the 
statistically significant features are seen in more than one detector or spectral 
order where available. We conclude that, to date, no credible X-ray line feature 
has been detected in a GRB afterglow. Finally, in a majority of cases, we find no 
evidence for significant absorbing columns local to the GRB host galaxy, implying 
there is little evidence from X-ray observations that GRB preferentially explode 
in high-density environments. 

Subject headings: gamma rays: bursts — gamma rays: observations — methods: 
data analysis — methods: statistical — X-rays: general 



1. Introduction 

Detections of emission and absorption features in GRB X-ray afterglows have been 
claimed in a number of observations. Iron emission lines from neutral as well as ionized 
material have been reported in the spectra of GRB970508 (Piro et al. 1999a), GRB970828 
(Yoshida et al. 1999, 2001), GRB991216 (Piro et al. 2000), and GRB000214 (Antonelh et al. 
2000). Reeves et al. (2002) claim detections of several emission lines from mid-Z elements 
in the afterglow spectrum of GRB011211, observed with XMM-Newton. Similar features 
have been reported in a selected time interval of the XMM-Newton spectrum of GRB030227 
(Watson et al. 2003a), and in the Chandra high-resolution spectrum of GRB020813 (But- 
ler et al. 2003b). Watson et al. (2002a) also argue that the XMM-Newton spectra of both 
GRB001025 and GRB010220 are better fit by a thermal emission line model in collisional 
ionization equilibrium compared to a smooth power-law continuum model. If these inter- 
pretations are correct, line identifications, measurements of their equivalent widths, velocity 
shifts, and their temporal behavior provide extremely useful diagnostics for inferring the 
physical conditions, including the geometry and dynamics of the burst environment, as well 
as the nature of the progenitor star (see, .e.g, Lazzati, Campana, & Ghisellini 1999; Rees & 
Meszaros 2000; Weth et al. 2000; Paerels et al. 2000; Meszaros & Rees 2001; Lazzati & Perna 
2002; Lazzati, Ramirez-Ruiz, & Rees 2002; Ballantyne & Ramirez-Ruiz 2001; Ghisellini et 
al. 2002; Kallman, Meszaros, & Rees 2003; Kumar & Narayan 2003; see, also, Boettcher 
2003 for a recent review). 
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The statistical significances and, hence, the reahty of these detections are subject to 
confirmation. All of the features observed to date were found by the original authors to be 
significant at the ~ Sa level, assuming that the line identifications and their energies are 
known a priori (i.e., single trial), with the exception of an emission line in the GRB991216 
spectrum, reported to be detected at 4.7cr (Piro et al. 2000). In the case of GRB011211, 
where the detection of multiple emission lines is claimed, the significance is estimated to 
be as high as 99.97 — 99.98% based on Monte Carlo simulations and Ax^ tests (Reeves 
et al. 2003, 2002). Borozdin & Trudolyubov (2003), however, criticize the data reduction 
and background subtraction procedures performed by Reeves et al. (2002), arguing that the 
single-trial significance of even the strongest feature is only S.lcr. Rutledge & Sako (2003) 
criticize the statistical analysis method adopted by both Reeves et al. (2003) and Reeves et 
al. (2002), concluding that the lines would not be discovered in a blind search. Interestingly, 
Watson et al. (2002b) report that the afterglow spectrum of GRB020322, one of the brightest 
afterglows observed with XMM-Newton, does not show any emission lines. The Chandra 
Low-Energy Transmission Grating (LETG) observation of GRB020405 (Mirabal, Paerels, 
& Halpern 2003), as well as the well-exposed High- Energy Transmission Grating (HETG) 
observation of GRB021004 (Sako & Harrison 2002a; Butler et al. 2003b) also appear not to 
show evidence for any strong discrete spectral features. 

Given the low significance of the individual detections, and the varied techniques used 
to analyze the spectra, it is difficult to assess the reality of these features from the literature 
alone. The primary problem is that the line energies are a priori unknown, and this is often 
not properly accounted for in deriving the true statistical significance. There is no consistent 
theoretical model which unambiguously predicts the line energies. In some cases, features 
are interpreted as near-neutral Fe Ka fluorescent emission, sometimes as highly-ionized Fe 
recombination lines, and others report thermal line emission from mid-Z elements with 
abundances often adjusted to super-solar values (in the case of thermal models) to match 
observations. In most cases, arbitrary blueshifts are invoked to adjust the closest atomic 
transition to match the observed energy of detected excesses (Reeves et al. 2003, 2002). In 
cases where the GRB redshift is not known, it is chosen to match line energies to probable 
atomic transitions. This means that the line energies are, in fact, free parameters, and this 
must be corrected for in deriving the detection significance. 

A secondary issue is the statistical method applied to argue for the detection of lines. 
A large fraction of the claims in the literature are based on using the F-test. This test 
compares the x^ values derived from two different models (e.g. with and without lines), and 
uses an analytic statistical model to determine if the difference represents an unusually good 
improvement. A probabilistic confidence determines if the improvement implies one model is 
more likely than the other to be correct. However, Protassov et al. (2002) point out that the 



-4- 



analytic F-test is not applicable to emission line detection, and will give both false positives 
and false negatives. Claims of line detections based on the F-test are therefore subject to 
scrutiny. 

The goal of this paper is to present a comprehensive, uniform data analysis of high- 
quality X-ray afterglow spectra. Numerous theoretical papers have been written on the 
interpretation and modeling of claimed emission line detections, and we do not attempt to 
validate or dispute specific interpretations. Rather, we perform our analysis independent 
of any particular model. We present detailed spectral analyses of 14 previously published 
datasets, in addition to 7 observations for which no detailed analysis has been published in 
a refereed journal. 

The paper is organized as follows. In §2, we present our sample and briefly summarize 
some of the relevant properties and previous observational analyses of each of the bursts. 
In §3, we describe the data reduction procedures we adopt for the various instrument con- 
figurations used in the observations. Results from the spectral fits are presented in §4. We 
describe our Monte Carlo simulations in §5. The results are discussed in §6 along with a 
brief comment on each of the sources with reported emission line detections. Finally, we 
discuss some of the implications of our results in §7. 



2. Observations 

The sample we have analyzed consists of 21 GRB X-ray afterglows up to and including 
GRB040106, and includes the brightest sources that allow meaningful constraints on discrete 
spectral features, and all sources for which the detection of emission lines has been reported. 
The source list and some of the relevant information are summarized in Table 1. Below we 
provide a brief description of the X-ray observations and previous analyses of each event, 
along with other outstanding features. 

GRB970508. Initially detected by the Beppo-SAXWFl (Costa et al. 1997), GRB970508 
is the second burst with an optical afterglow detection, and was the first event with a 
measured redshift {z = 0.835; Bloom et al. 1998). The Beppo-SAX NFl observed the error box 
starting UT May 9.1375 (0.43 days after the burst), detecting a relatively bright {F2-10 kcv ~ 
10"^^ erg cm^^ s^^) X-ray afterglow. The Beppo-SAX NFI continued to monitor the X-ray 
afterglow, performing a total of four observations, the last taking place six days after the 
burst (Piro et al. 1999a). 

Piro et al. (1999a) reported a narrow line feature at 3.5 keV, corresponding to Fe XXVI 
Lya line emission at the measured host redshift, during part of the afterglow. The afterglow 
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did not decay smoothly, but the X-ray hghtcurve as measured by the NFI shows a bump, 
beginning a day after the event, lasting about two days. By breaking the first observation into 
two parts (denoted "la" and "lb"), prior to and after the "bump", Piro et al. (1999a) find 
evidence for line emission in the first segment at 99.3% (2.7o"; using an F-test) confidence, 
assuming the known redshift. No line emission is found in the second part of the pointing, 
or in any of the subsequent observations. 

GRB970828. This event is notable for being the prototype dark GRB. A fading X-ray 
afterglow was discovered, but no associated optical transient was found despite extensive 
followup. Initially positioned by the RXTE AW Sky Monitor (ASM; Remillard et al. 1997), 
ASCA followed up the location beginning 1.2 days after the event, detecting a moderately 
bright {F2-10 kcv = 4 X 10^^^ erg cm"^ s"^) X-ray afterglow (Yoshida et al. 1999). Like 
GRB970508, GRB970828 exhibited variability in the X-ray superimposed on the power-law 
decay. 

Dividing the observation into three segments, denoted A, B, and C, Yoshida et al. 
(1999) found an excess above a power-law spectrum centered at 5.04to3i keV with a width 
of 0.31^Q3^ keV in segment B with 98.3% confidence (2.4o"; using an F-test), which coincided 
with flaring activity in the afterglow. Both A and C have spectra consistent with an absorbed 
power law. Initially Yoshida et al. (1999) identified the line with Fe Ka fluorescence emission, 
however the subsequent measurement of the redshift of 0.9578 for the host (Djorgovski et al. 
2001) places the rest-frame energy at 9.87 keV, forcing Yoshida et al. (2001) to re-identify 
the feature with an Fe XXVI recombination continuum with no associated line emission. 

GRB991216. This event was a very bright BATSE trigger (Kippen, Preece, & Giblin 
1999). The PGA on RXTE scanned the error circle, and discovered a bright (5 mCrab) 
X-ray source identified as the GRB afterglow (Takeshima et al. 1999). The discovery of an 
associated optical transient (Uglesich et al. 1999) lead to a redshift measurement {z = 1.02; 
Vreeswijk et al. 1999). Chandra followed up the optical afterglow position with the HETG 
in a 9.65 ksec observation 37 hours after the burst, detecting a bright (-F2_io keV = 2 x 10^^^ 
erg cm~^ s^^) X-ray afterglow (Piro et al. 1999b). 

Piro et al. (2000) analyzed the HETG spectrum, and found excess emission at two 
energies, which they interpret as due to highly ionized iron. The most prominent feature, 
which they claim to be significant at the 4.7a level, lies at 3.49 ± 0.06 keV and is much 
broader than the instrument resolution, with a gaussian width of o" = 0.23 ± 0.07 keV. This 
is, by far, the highest significance claimed for any emission line of any GRB. The second, 
marginally significant feature lies at 4.4 ± 0.5 keV (99.5% confidence, or 2.8o"; F-test). If 
associated with the iron recombination continuum (threshold rest energy at 9.28 keV), the 
implied redshift isz = l.ll±0.11, and if the lower energy feature is associated with H-like 
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Lya of Fe, the implied z = 1.0 ± 0.02. These are both consistent with the measured z of 
1.02. 

GRB000210. This is another event categorized as a dark GRB. The Beppo-SAX WFC 
positioned this burst, which was among the brightest detected by that mission (Stornelh 
et ah 2000). Deep optical searches performed ~ 16 hours after the event failed to locate 
an optical afterglow (Gorosabel et al. 2000b). Both the Beppo-SAX NFl as well as Chandra 
ACIS-S performed followup observations, and both detected the X-ray afterglow; Beppo-SAX 
observations started ~8 hours after the event (Costa et al. 2000), and a Chandra pointing 
began 21 hours after (Garcia et al. 2000). Deep optical observations of the C/ian(/ra position 
revealed a host galaxy, with a measured redshift of 0.8463 (Piro et al. 2002). A radio afterglow 
was also detected (Piro et al. 2002), and a ~ 2.5a sub-mm excess toward the GRB000210 
has been reported (Berger et al. 2003a). 

The X-ray afterglow appears unremarkable, following a smooth t^^-^^ decay. Fits to the 
joint 'HYl/ Chandra data reveal a spectrum consistent with an absorbed power law (Piro et 
al. 2002). 

GRB000214. This burst was detected by the BeppoSAX G^YiM and positioned by the 
WFC (Piro 2000a,b). Followup with the NFI beginning 12 hours after the burst revealed 
an X-ray afterglow with a flux of -^2-10 kcv ~ 8 x 10~^^ erg cm~^ s~^ (Antonelli et al. 
2000). The X-ray lightcurve is unremarkable, following a power-law decay. No optical 
counterpart was found, however the followup was complicated by the moon, and the limits 
are not particularly constraining. No redshift has been derived due to the lack of an accurate 
localization, required for host identification. 

Antonelli et al. (2000) fit the X-ray spectrum of the entire NFI observation. They find 
an absorbed power law to be a poor fit (chance probability P = 0.02), with an excess above 
the power law at an energy between 4 — 5 keV. By adding a narrow gaussian line, they obtain 
an acceptable fit, with the line centroid being at 4.7 ± 0.2 keV, with an equivalent width of 
~ 2.1 keV, with 99.73% confidence {3a, using an F-test). Although no independent redshift 
was measured Antonelli et al. (2000) interpreted this line as a redshifted Fe Ka line. 

GRB000926. The Interplanetary Network (IPN) discovered GRB000926 (Hurley et al. 
2000). The optical afterglow of this 25 s long event was identified less than a day later 
(Gorosabel et al. 2000a; Dall et al. 2000). The redshift, measured from optical absorption 
features, is 2.0379 ± 0.0008 (Castro et al. 2003). The afterglow was well-monitored in the 
optical (Harrison et al. 2001b), and was detected in the IR (di Paola et al. 2000; Fynbo et 
al. 2001) and radio (Harrison et al. 2001b). 

The Beppo-SAX Medium Energy Concentrating Spectrometer (MFCS) instrument (sen- 
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sitive from 1.3 — 10 keV) discovered the X-ray afterglow (Piro et al. 2001). The X-ray source 
was weak, and so we do not analyze the Beppo-SAX data. Garmire et al. (2000) observed 
the source for 10 ksec with Chandra using the ACIS-S3 backside-illuminated chip on Sept 
29.674 - 29.851. Chandra again observed GRB000926 in a 33 ksec long ToO taken Oct 
10.176 — 10.76 UT, also with ACIS-S3. The afterglow was clearly detected in each of these 
observations. Piro et al. (2001) fit the combined Chandra and NFI spectrum, finding an 
absorbed power law to be an acceptable fit, with no evidence for any spectral features. 

GRB001025. This burst (also called GRB001025A, to differentiate from a second GRB 
observed during the same UT day), which lasted about 15 s, was localized by the RXTE 
ASM (Smith et al. 2000). No optical counterpart was identified to a limit oi R = 24.5 
(Fynbo et al. 2000). XMM-Newton observed the location of GRB001005 starting about 1.9 
days after the burst with a total EPIC exposure of about 25 ksec, and found two X-ray 
sources in the ASM error circle (Altieri et al. 2000). The brighter of the two, detected at 
-^0.2-10 kcv = 4.4 X lO"^"' erg cm~^ s~^, is considered the most likely afterglow candidate, 
however some uncertainty remains over this identification. Watson et al. (2002a) find the 
source to have decreased in flux during the observation at 99.8% confldence, while Borozdin 
& Trudolyubov (2003) flnd the source to be consistent with constant flux. The lack of 
an optical counterpart to the X-ray source to R = 24.5 would tend to rule out an AGN, 
strengthening the association with the GRB. 

Borozdin & Trudolyubov (2003) analyzed the spectrum of the afterglow candidate, and 
found an acceptable flt using a power law with Galactic absorption. Watson et al. (2002a), 
however, flnd a collisionally ionized plasma model to be a signiflcantly better flt to the 
spectrum. In the thermal plasma model, the redshift is allowed to vary to obtain the best 
flt. According to Watson et al. (2002a), the plasma model is a better flt due to discrete 
fluctuations observed between ~ 0.5 and 2 keV, which these authors interpret as due to the 
blend of lines from Mg, Si, S, Ar, and Ni L. The flt temperature of 3.4 keV is determined 
by the continuum shape, and the redshift corresponding to the best flt is 2; = 0.70]to'i^, 
from reported observed line energies (with detection signiflcances of the individual lines) of 
O.SOtom (4 X 10"^), 1.16±0.05 (8 x 10"^), 1.64±0.07 (8 x 10"^), 2.2±0.1 (0.02), and 4.7+^;^ 
keV (0.12), corresponding to Mg XII (1.46 keV), Si XIV (1.99 keV), S XVI (2.60 keV), Ar XVII 
(3.30 keV), and Ni XXVIII (8.10 keV), respectively. While the probabilities for detecting the 
individual lines were given, the method for deriving these probabilities was not. 



^Watson et al. (2002a) give this redshift in the text; but in their Table 1, they state the best-fit is z = 0.53 
(90% confidence range 0.5-0.55). Our own fit to the given line centroids with uncertainties using the line 
identifications of Watson et al. (2002a), which does not include contimuim uncertainties, gives z — 0.65±0.04. 



GRB010220. This burst, positioned by Beppo-SAX (Piro 2001a,b; Manzo et al. 2001), 
has no identified optical counterpart. Searches reached to 23.5 in R, but were comphcated by 
the position of the burst in the Galactic plane (Berger et al. 2001). XMM-Newton observed 
the error circle starting 14.8 hours after the event. Like GRB001025, the X-ray counterpart 
is uncertain, with four sources discovered in the error circle. None of these is clearly variable, 
and the brightest has been associated with the afterglow (Watson et al. 2002a). 

Watson et al. (2002a) have analyzed the spectrum of the probable counterpart, and 
find an absorbed power law to be an acceptable fit. They also claim a significant deviation 
from a power-law model at around 3.9 keV, and that adding an unresolved gaussian line 
with equivalent width l.S^Ji ^^^ improves the fit at > 99% {2.6a) confidence (as with 
GRB0001025, this reference gives the confidence level for detection of this line, but does not 
describe the method by which the confidence level was determined). According to their fits, 
a thermal plasma model is significantly preferred over the absorbed power law, and allowing 
the Ni abundance to vary from solar, the 3.9 keV can be accounted for with this element. 
It is unclear how Watson et al. (2002a) identify Ni (rest energy 7.47 keV) with the feature, 
since no optical redshift has been measured. The required redshift assuming the 3.9 keV line 
to be due to Ni Ka fiuorescence is z = 0.92. 

GRB011030. This X-ray rich GRB (or X-ray fiash - XRF) was originally detected by 
Beppo-SAX (Gandolfi 2001a,b). While it is not clear whether XRFs and GRBs are different 
manifestations of the same phenomenon, XRFs exhibit afterglows with similar characteristics 
to GRB afterglows. A radio afterglow was detected ~ 9 days after the burst (Taylor, Frail, 
& Kulkarni 2001). No optical afterglow was detected for this event. Chandra performed 
two followup observations, 10.5 (Harrison et al. 2001a), and 20 days after the event. With 
ACIS-S, Chandra detected an X-ray source with a fiux of -F2-10 kev2.4 x 10^^^ erg cm^^ s^^ in 
the error box, which faded by more than an order of magnitude between the two pointings. 
This source is certainly the afterglow associated with the initial fiash. The spectrum, on 
preliminary analysis, appears consistent with a power law with Galactic absorption (Harrison 
et al. 2001a). 

GRB011211. This X-ray rich GRB was localized by Beppo-SAX, and is notable, at 270 s 
duration, for being the longest event detected by the WFC (Gandolfi & Halpern 2001). The 
identification of an optical afterglow lead to the measurement of the absorption-line redshift 
oi z = 2.140 (Holland et al. 2002). XMM-Newton began observing the error circle starting 
~11 hours after the event, and discovered an X-ray afterglow (Santos-Lleo et al. 2001) with 
a time-averaged fiux of -Fq. 22-10 keV = 1-7 x 10^^^ erg cm~^ s^^ (Reeves et al. 2003). 

Considerable controversy surrounds the analysis of the X-ray spectrum of this event. 
Reeves et al. (2002) consider only the PN data, divide the observations into time sections. 
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and claim that in the first 5 ksec interval, an absorbed power law marginally fits the data. 
They find that adding lines from partially ionized Mg, Si, S, Ar and Ca, with a mean redshift 
of 1.9 (implying a blueshift, possibly due to an outfiow, of ~ 0.1c) improves the fit at 99.5% 
confidence. The corresponding observed (rest-frame) line centroids are 0.44, 0.71, 0.88, 1.22 
and 1.46 keV. They find the most significant of the lines (0.71, 0.88, and 1.22 keV) are detected 
with 98% (2.3(t), 97% (2.2o"), and 93% (1.8cr) confidence for detection at a single energy, 
respectively, employing the F-test. Borozdin & Trudolyubov (2003) analyze the same 5 ksec 
interval both for the PN data alone, and for the combined PN, MOSl and M0S2 data. For the 
combined set they find a good fit to an absorbed power law with Galactic absorption, with no 
improvement from adding lines. For the PN data alone, they find marginal improvement for 
adding lines at the energies found by Reeves et al. (2002), and they suggest the difference is 
due to the selection of background regions. Rutledge & Sako (2003) independently analyzed 
the statistical significance of the features, finding that significance is marginal. 

GRB020321. This weak GRB was discovered in the Beppo-SAX WFC using ground 
trigger logic (Guidorzi et al. 2002). A very faint variable object was detected in a foUowup 
observation with the NFI ~6 hours after the event (Gandolfi 2002a). Deep optical imaging 
failed to find any variable source either at the NFI position, or the WFC error box to a 
limiting magnitude of i? ~ 24 (Price, Dressier, & McCarthy 2002). A refined analysis of 
the NFI data showed the detection to be marginal, and an XMM-Newton observation taking 
place 10.3 — 24.2 hours after the GRB failed to detect the NFI source, but discovered another 
variable X-ray source in the WFC error circle. This source was not detected in a Chandra 
observation taking place 9.9 days after the GRB, and is considered the likely (although not 
certain) counterpart of the GRB (In't Zand et al. 2002). We analyze the data from this 
source in this paper. 

GRB020322. This weak Beppo-SAX hnrst was followed up with the NFI 6.5 hours after 
the event, and a new, relatively bright X-ray source was found in the error circle (Gandolfi 
2002b). This was confirmed, in XMM-Newton observations starting about 15.5 hours after 
the GRB to be the fading afterglow. The X-ray flux, at -Fo.2-10 kcv = 3.5 x 10"^^ erg cm~^ s~^ 
made this a high signal-to-noise XMM-Newton detection. A very weak optical counterpart 
was detected (Bloom et al. 2002), and confirmed to be the afterglow. The host, at 27th 
magnitude (Burud et al. 2002) is very faint, and no redshift has been measured. Watson et 
al. (2002b) analyzed the spectrum, and find it to be consistent with an absorbed power law, 
with column in excess of the Galactic value. 

GRB020405. This bright IPN burst (Hurley et al. 2002) was very well studied, with 
bright optical (Price et al. 2003a) and radio (Berger et al. 2003b) transients, evidence for 
a late-time red bump in the optical decay possibly due to a supernova component Price et 
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al. (2003a), optical polarization detections (Covino et al. 2003), and a measured redsliift of 
z = 0.690 (Masetti et al. 2002; Price et al. 2003b). Detailed optical/NIR spectrophotometric 
and polarimetric analyses were published by Masetti et al. (2003). Mirabal, Paerels, & 
Halpern (2003) triggered a long 50 ksec Chandra LETG observation, which started 1.68 
days after the GRB. They detected an X-ray afterglow with flux Fo.2-10 kev = 1-36 x 10~^^ 
erg cm^^ s^^ 1.71 days after the event. They analyze the X-ray spectrum, finding it to be 
consistent with a featureless power law with the suggestion of an absorbing column in excess 
of Galactic. 

GRB020813. This typical long-duration GRB was localized by HETE-2 (Villasenor et 
al. 2002). The rapid dissemination of the position enabled prompt identification of an optical 
counterpart (Fox, Blake, & Price 2002), and the redshift z = 1.255 has been measured both 
in absorption and in emission (Barth et al. 2003). Chandra observed the region with the 
HETG starting 17.7 hours after the burst, with a total integration of 77 ksec (Vanderspek 
et al. 2002). The mean flux of the X-ray afterglow during this observation was -Fo.6-6 kcv = 
2.2 X 10^^^ erg cm^^ s^^. A preliminary reduction showed no obvious spectral features 
(Vanderspek et al. 2002). Subsequently, Butler et al. (2003b) reported evidence at 3.3o" 
significance for a line at 1.3 keV, which they interpret as due to H-like sulfur, blueshifted 
(with respect to the host redshift) by 0.1c. 

GRB021004. The position of this HETE-2 GRB was disseminated quickly (Shirasaki 
et al. 2002), and a bright optical afterglow discovered 9 minutes after the event (Fox et al. 
2003b). The redshift was soon determined to be z = 2.323 from the identification of the Lya 
line (Chornock & Filippenko 2002). The Chandra HETG observed the afterglow starting 20.5 
hours after, with a total exposure of 87 ksec. Preliminary reduction of the data showed a 
fading X-ray source with mean flux of -F2-10 kcV = 4.3 x 10^^^ erg cm^^ s^^ (Sako & Harrison 
2002b). This analysis also showed the spectrum to be consistent with a power law, with no 
evidence for absorption in excess of the Galactic value. 

GRB030226. An optical afterglow of this average brightness, 100 s long HETE-2 burst 
(Suzuki et al. 2003) was discovered 2.6 hours after the trigger (Fox, Chen, & Price 2003a). 
Absorption spectroscopy identified the redshift as z = 1.99 (Ando et al. 2003). 

Chandra observed this event with ACIS-S3 starting 37.1 hours after the trigger, with a 
total exposure of 34 ksec. A relatively weak; ^2^10 kcV = 3.2 x 10^^^ erg cm^^ s^^ averaged 
over the observation. X-ray afterglow was clearly detected at the position of the optical 
transient. A preliminary analysis of the spectrum (Pedersen et al. 2003) found it to be 
consistent with an absorbed power law, with absorption in excess of Galactic. A subsequent 
analysis disputed the Nh value, finding it consistent with the Galactic value (Sako & Fox 
2003). 
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GRB030227. This weak, long-duration GRB was positioned by Integral IBIS (Gotz, 
Borkowski, & Mereghetti 2003a; Gotz, Mereghetti, & Borkowski 2003c), and followed up 
in the X-ray by XMM-Newton with a 13 hour pointing beginning 8 hours after the event 
(Mereghetti et al. 2003). After the detection of a fading X-ray source with average flux of 
-fo.2-10 kcv = 8.7 X 10"^^ erg cm~^ s~\ a relatively faint optical counterpart was discovered in 
the XMM-Newton error circle (Soderberg et al. 2003a). No redshift has yet been measured. 

Mereghetti et al. (2003) analyzed the spectrum from the entire observation, and find it 
to be consistent with an absorbed power law, with a column in excess of the Galactic value. 
They also find weak evidence for an emission feature at 1.67 keV, which they postulate could 
be due to iron, without asserting a detection. In a subsequent analysis, Watson et al. (2003a) 
divided the observation up into four segments. In the first three, the spectrum in consistent 
with an absorbed power law, but in the last 11 ksec, they claim the presence of 5 emission 
lines at centroid energies of 0.62, 0.86, 1.11, 1.35, and 1.66 keV, with significance ranging 
from 1.7 — 3.8(7, using an F-test, for detection of lines at single energies depending on how 
they estimate the uncertainty. Assuming a redshift of 1.35, the line energies correspond 
to H- and He-like lines of Mg, Si, S, Ar, and Ca. Watson et al. (2003a) also claim that 
the continuum fit for zero emission-line fiux is ruled out at the 4.4 — 4.7(T significance. For 
comparative purposes, Mereghetti et al. (2003) estimate the significance of the feature at 

0.01 
0.03 



E = 1.67+o:o3 keV to be 3.2a using an F-test 



GRB030328. The position of this typical long-duration HETE-2 GRB was rapidly 
disseminated (Villasenor et al. 2003), enabling both the early identification of an optical 
transient, and a Chandra grating observation beginning only 15.3 hours after the initial 
trigger (Peterson & Price 2003). Optical spectroscopy identified the source redshift as 2; = 
1.52 (Martini, Garnavich, & Stanek 2003). 

The LETG observed the event for 94 ksec. A preliminary analysis (Butler et al. 2003a) 
found that during the observation the mean source fiux was -Fo.5-3 keV = 1.9x 10~^^ erg cm~^ s^ 
and that the spectrum is consistent with a power law with Galactic absorption. 

GRB030329. The rapid dissemination of the position of this HETE-2 event (Van- 
derspek et al. 2003) enabled rapid identification of a bright optical counterpart, and the 
measurement of the redshift, z = 0.169 (Greiner et al. 2003). The proximity and brightness 
of this famous burst enabled the detection of a Type Ic supernova, SN2003dh, superimposed 
on the afterglow, securing the association of some long-duration events with the deaths of 
massive stars (Garnavich et al. 2003; Stanek et al. 2003; Hjorth et al. 2003). The X-ray 
afterglow was bright enough to be detected by RXTE approximately 5 hours after the burst 
(Marshall & Swank 2003). XMM-Newton observed the position twice, in pointings separated 
by nine days. Temporal and spectral analyses of the first observation and an analysis of the 
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temporal behavior through the second observation were pubhshed by Tiengo et al. (2003) 
and Tiengo et al. (2004), respectively. The first observation had a 30 ksec exposure, while 
the second was 39 ksec. No detailed line search analysis has yet been published of these 
data. 

GRB031203. This long burst was localized by IntegrallBlS (Gotz et al. 2003b), and 
a followup observation by XMM-Newton started only ~ 6 hours after the burst. Of the 
two bright sources that were detected inside the IBIS error circle (Rodriguez-Pascual et al. 
2003), one of them appears to be not consistent with the RO SAT upper hmit, suggesting 
that this is the afterglow (Campana et al. 2003). The redshift of the possible host galaxy 
was determined to be z = 0.105 (Prochaska et al. 2003), coincident with a variable radio 
source (Frail 2003; Soderberg, Kulkarni, & Frail 2003b). A spectacular X-ray dust echo was 
detected for the first time in a GRB (Vaughan et al. 2003, 2004). 

A second XMM-Newton observation was performed starting approximately ~ 2.9 days 
after the burst. Watson et al. (2004) claim that the X-ray light curve steepens from — 0.4±0.1 
to —0.9 ± 0.1 approximately one day after the burst, although we cannot confirm this in our 
independent analysis. Both observations appear to be consistent with a single power-law 
decay time slope of —0.4 ± 0.1, which is one of the most gradual decay ever detected. 

GRB040106. XMM-Newton observed this long GRB localized by IntegrallBlS (Mereghetti 
et al. 2004) starting ~ 6 hours after the burst. A fading X-ray source was soon discovered 
(Ehle, Gonzalez- Riestra, & Gonzalez- Garcia 2004). This burst was the brightest ever ob- 
served by XMM-Newton, with ~ 23000 total photons detected in the 0.2 — 10 keV range. 
A likely optical counterpart has been identified (Masetti et al. 2004), but no spectroscopic 
redshift has been measured to date. 

A summary of reported emission line detections, including their identifications and the 
quoted significances, is listed in Table 2. This table lists two additional sources (GRB001025A 
and GRB010220) that are not listed in Table 1 of Boettcher (2003). 



3. Data Reduction 

In this section, we describe the reduction procedures that we adopted for the various 
observatory/instrument configurations. Table 3 lists, for each dataset, the total number 
of counts detected inside the extraction region (source plus background) and the estimated 
number of background counts. In a few cases, we split the observation into two time intervals 
to match those adopted in previous publications. 
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3.1. Chandra 

Three sources (GRB991216, GRB020813, and GRB021004) were observed with the Chan- 
dra (Weisskopf et al 2000) HETG (Canizares et ah 2000) with the ACIS-S (Garmire et al. 
2003) placed at the focal plane. The HETG consists of two separate grating arrays; the 
Medium Energy Gratings (MEG) and the High Energy Gratings (HEG), which are opti- 
mized in the approximate wavelength ranges A = 5 — 20 A and A = 1.5 — 15 A, respectively. 
The spectral resolving power of the HEG is R = \/A\fwhm ~ 1000 x (A/10 A) and is approx- 
imately twice that of the MEG, for a total of ~ 1350 and ~ 700 useful resolution elements 
for the HEG and MEG, respectively. Two other sources (GRB020405 and GRB030328) were 
observed with the LETG (Brinkman et al. 2000) with the ACIS-S at the focal plane. This 
configuration optimizes the sensitivity in the A = 5 — 30 A region with a resolving power 
that is approximately less than half that of the MEG, with ^ 400 resolution elements. 

The data were retrieved from the Chandra Data Archive^ and were processed using 
CIAO v2.3^. We use sky coordinates of the aspect-corrected Level 1 events to determine 
the location of the peak X-ray flux in the zeroth order image, which we use to measure the 
dispersion angles and assign wavelengths to each of the dispersed events. Only events with 
grade=0,2-6 and status=0 were used in the analyses. Background events produced during 
several serial readout frames on ACIS-S4 were removed using destreak, which significantly 
reduces the noise level on the positive dispersed orders. Source events were then spatially 
extracted using a 4.8" filter in the cross-dispersion direction. We then use a standard pulse- 
height filter to separate the first order events. Ancillary response files were generated based 
on these extraction regions. The relative effective area may be uncertain by ~ 10% above 
~ 2 A and ~ 20% below ~ 2 A^ 

We use also the zeroth order spectrum, which has an intrinsic spectral resolution that 
is lower than that of the HETG by an order of magnitude or more, but contains roughly 
an equal number of photons compared to the combined first order dispersed events. Source 
photons are extracted from a circular region with a radius of 4". Background events are 
extracted from several large off-axis regions on the ACIS-S chips. The positive and negative 
first order events are combined, which are then fit simultaneously with the zeroth order 
spectrum. The relative normalizations of the zeroth and first order spectra are allowed to 
vary within 20%. 



^http://asc. harvard.edu/cda/ 
•^ http://asc.harvard.edu/ciao/ 
''http://asc.harvard.edu/proposer/POG/html/ 
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Four sources (GRB000210, GRB000926, GRB011030, and GRB030226) were observed 
with the ACIS-S at the focal plane with no gratings. The spectral reduction procedure 
is identical to that of the zeroth order image in the HETG/ ACIS-S configuration described 
above. The effective area is ~ 15 and ~ 7 times higher at ii^ = 1.5 keV and 6 keV, respectively, 
than that with the HETG placed in front. The spectral resolving power of ACIS-S is i? ~ 
-Ekev/(0.110 -|- 0.012 -Ekev), where i^kev is the energy in keV, for ^ 190 resolution elements in 
the 0.5 — 8 keV range. None of these four datasets has been published to date. 



3.2. X MM- Newton 

Nine afterglows (GRB001025, GRB010220, GRB011211, GRB020321, GRB020322, GRB030227, 
GRB030329, GRB031203, and GRB040106) were observed with XMM-Newton Observatory 
(Jansen et al. 2001). XMM-Newton consists of three X-ray telescopes with the European 
Photon Imaging Camera (EPIC) PN (Striider et al. 2001) at one of the focal planes and the 
MOS (Turner et al. 2001) CCDs at the remaining two. Behind the two mirror assemblies 
that focus light onto the MOS are two identical sets of Reflection Grating Spectrometers 
(RGS; den Herder et al. 2001). The EPIC PN and MOS consist of CCD arrays with spectral 

1/2 

resolving powers of ~ 15 -E^cV ^^^ both the PN and MOS for ^ 220 resolution elements in 
the 0.2 — 8 keV range. The PN chips are sensitive in the ~ 0.2 — 10 keV range with a peak 
effective area of ~ 1200 cm^ at ~ 2 keV. The MOS chips are sensitive in the ~ 0.5 — 10 keV 
with roughly half the effective area per array of the PN. Data obtained with the Optical 
Monitor (OM) are not used for any of our analyses. 

The data were retrieved from the XMM-Newton Science Archive^ and were processed 
with the XMM-Newton Science Analysis Software^ (SAS). We use only data collected on the 
PN and MOS detectors. The RGS data are not used, since the signal-to-noise ratios (S/N) of 
each of the bursts are too low to provide useful constraints. Events on bad/hot pixels and 
pixels on the edge of the CCD chips were excluded from the analysis by selecting events with 
only FLAG=0. Event grade selections were made using PATTERN<=4 and PATTERN<=12 for the 
PN and MOS, respectively. Source events from both the PN and MOS were extracted from 
circular region centered on the source with a radius of 30 — 45" depending on the level of the 
background flux during the observation. Background events were taken from large source- 
free regions (typically ~ 2 — 3' in radius) that lie on the same CCD chip. Redistribution 
matrices and effective area curves were computed for each observation and each instrument 



'''http://xmiii.vilspa.esa.es/external/xmni_data_acc/xsa/index.shtnil 
^http://xmni.vilspa.esa.es/external/xmni^w_cal/sas_franic.shtnil 
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using the SAS tools rmf gen and arf gen, respectively. The PN and MOS spectra are fit 
simultaneously, allowing the relative normalizations to vary within 5% 



v'o\ 



For GRB010220, Watson et al. (2002a) found 20 ksec (of 43 ksec exposure) in which the 
background flaring was sufficiently low for analysis, while find find only ~ 12 ksec. The data 
of GRB011211 are split into two time intervals; (1) the first 5 ksec during which detections 
of soft X-ray lines were reported by Reeves et al. (2002) and (2) the remaining 24 ksec where 
no lines were seen. The data of GRB030227 were analyzed using data collected over the 
entire observation, and using data only collected during the last 11 ksec, during the period 
reported by Watson et al. (2003a) to contain emission lines. GRB030329 was observed at two 
epochs separated by ~ 9 days, and the spectra are analyzed separately. For the remaining six 
sources (GRB001025, GRB010220, GRB020321, GRB020322, GRB0301203, and GRB040106), 
we use data collected over the entire observation. 



3.3. Beppo-SAX 

Two sources (GRB970508 and GRB000214) were observed with the Beppo-SAX NFI, 
which consists of the Low Energy Concentrator Spectrometer (LEGS; Parmar et al. 1997) 
and two working units of the Medium Energy Concentrator Spectrometer (MEGS; Boella 
et al. 1997). The LEGS is sensitive to X-rays in the range 0.1 — 4.0 keV with an effective 
area of 40 cm^ at i? ~ 2 keV. The MEGS units are sensitive a slightly higher energy range 
(1.3 — 10 keV) with a combined peak effective area of 90 cm^ at -E ~ 5 keV. The spectral 
resolving power of both the LEGS and the MEGS can be represented approximately by 
i? = 10 (-E/6.4 keV)°'^^, for ^ 15 energy independent resolution elements (0.1 — 3.0 keV, 
LEGS; Parmar et al. 1997) and ~ 60 resolution elements (1.3 — 8 keV, MEGS; Boella et al. 
1997). 

The data reduction was performed with HEASOFT v5.2^. The MEGS and LEGS events 
were extracted from circular regions with radii of 3' and 8', respectively, centered on the 
source. Since the background depends on the detector location, background spectra were 
generated from deep exposures of blank fields using an identical extraction region. The 
background flux levels were normalized by comparing the count rates in the outer source- 
free regions. For both observations, the MEGS background fluxes were found to be within 



^see, e.g., http://xmm.vilspa.esa.es/docs/docunients/CAL-TN023.ps.gz 
^http://hcasarc. gsfc.nasa.gov/docs/software/lheasoft/ 
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5%. Pre-made response matrices and effective area files^ were used for the analyses. 

GRB970508 was observed with Beppo-SAX during muhiple epochs within one week of 
the burst (Piro et al. 1998). The presence of an intense iron hne was reported during its early 
afterglow (< 1 day; time interval denoted as "la" in Piro et al. (1999a)), and disappeared 
after ~ 1 day (interval denoted as "lb"). We analyze these spectra separately. 



3.4. ASCA 

The afterglow of GRB970828 was observed by the ASCA Observatory (Tanaka, Inoue, 
& Holt 1994) and the data were originally published by Yoshida et al. (1999). ASCA carries 
four focusing mirrors, with two CCD imaging spectrometers (SISO and SISl) and two gas 
scintillation imaging spectrometers (GIS2 and GIS3) at the focal plane. The data were 
retrieved from the HEASARC archive^'^. 

Yoshida et al. (1999) identify a feature during a ~ 20 ksec interval of flaring activity, 
which is denoted as time segment "B" by Yoshida et al. (1999). The line was claimed to 
be absent during the other time intervals. We, therefore, sum the remaining data (time 
segments "A" and "C") and analyze them collectively. 

The data reduction was again performed using HEASOFT v5.2. The data were screened 
through standard criteria: rejection of hot flickering pixels and data contaminated by bright 
earth, grade selections, etc. Source events were extracted from circular regions with a radii of 
3' and 5' for the SIS and GIS, respectively. The background flux level was steady during the 
entire observation so we use the total exposure to create background spectra for each instru- 
ment using off-axis source-free regions. The data from all four units were flt simultaneously 
in the energy ranges 0.5 — 10 keV and 0.7 — 10 keV for the SIS and GIS, respectively. 



4. Spectral Analysis 

Our primary method for searching for line features is through Monte Carlo simulations, 
however for general spectral characterization, and for comparison with other work, we have 
performed standard spectral analysis, fltting the raw spectrum to an absorbed power law 
model. 



^ftp://ftp. asdc.asi.it/pub/sax/cal/responses/ 
^°http://heasarc. gsfc.nasa.gov/db-perl/W3Browse/w3browse. pi 
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We performed all of the spectral fitting with the XSPEC v.ll^^ spectral fitting package. 
We binned the spectra for all sources so that each spectral bin contains at least 15 - 20 
counts (the number depending on the total counts detected) except for the grating spectra, 
which were binned uniformly by a certain factor so that each bin typically contains at least 
a few counts on average^^. We then fit each of the background subtracted spectra using 
a power-law model with both Galactic and intrinsic absorption by cold material. For the 
GRB with known redshifts, we fix the intrinsic absorber to be at the measured redshift. 
Otherwise, we set z = 1, which roughly represents an average value of all GRB with known 
redshifts. The model contains three free parameters; the intrinsic column density N^^^, 
the power-law photon index F, and its normalization -Fikev defined as the flux in units of 
photons cm~^ s~^ keV~^ at 1 keV. The functional form of the continuum model can be 
written as, 

F{E) = Fikev i?k"ev e-^"^""" -ab.(i^,.)+^°^' ^.UE))^ (1) 

where £^kev is the photon energy in units of keV and crabs(-E) is the energy-dependent ab- 
sorption cross section of interstellar material. The Galactic hydrogen column densities are 
derived from H I maps by Dickey & Lockman (1990). We assume that the metal abundances 
relative to hydrogen are those given by Anders & Grevesse (1989). 

We manually incorporate a fix to the low-energy quantum efficiency degradation of the 
Chandra ACIS-S using the acisabs model^^. This effect, due to accumulation of molecular 
contamination onto either the optical filters and/or the CCD chips, is time- dependent and 
reduced the E ~ 0.7 keV effective area by ~ 20-30% during the GRB020813 and GRB021004 
observations and by ~ 40% during the GRB030328 observation compared to pre-launch 
values. The GRB991216 observation was performed only ~ 5 months after launch and the 
data are therefore not expected to be affected significantly. 

Figures 1-32 show the binned spectra, and x value for each spectral bin for the 
respective best-fit model. Table 3 lists the best-fit parameters and the resulting x^ ^''^^ 
null-hypothesis probability, P{x, ^)- We do not list P{x, ^) for the dispersed grating spectra 
because they were fit with a C-statistic due to the low number counts detected per spectral 
bin. The sample as a whole is consistent with the hypothesis of absorbed power law X-ray 



^^http://hcasarc. gsfc.nasa.gov/docs/software/lheasoft/xanadu/xspec/indcx. html 

^^We choose not to bin the grating spectra to contain at least ^ 15 counts per bin, since this typically 
results in an undersanipling of the instrument resolution by an order of magnitude or more in most parts of 
the spectra. We note, however, that the continuum parameters derived from spectra of both binning schemes 
are nearly identical and any subsequent analyses and their results are not affected by our choice of binning. 
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http://www.astro.psu.edu/users/chartas/xcontdir/xcont.html 
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emission; out of 27 separate spectra that are fit using x^-minimization, four have P{x, ^) < 
0.1, consistent with what one would expect statistically from a uniform sample. If we consider 
the four spectra with P{x, v) < 0.1; the zeroth order spectrum of GRB991216, GRB000214, 
the zeroth order spectrum of GRB021004, and the zeroth order spectrum of GRB030328 
(Figures 6, 8, 22, and 27, respectively), only GRB991216 and GRB000214 appear to show 
discrete residuals suggestive of statistically significant line features. Examination of the 
other spectra shows one other interesting case, GRB970828 segment B (Figure 4), where the 
residuals appear to be concentrated in a broad feature (although P{x, J^) is acceptable). As 
we discuss in the following section, however, the x^ statistic is not sensitive to weak and/or 
localized fluctuations, so the x^ values are not in general particularly useful for line searches. 
Similarly, F-tests yield unreliable values for the significances of line-like features as described 
in detail by Protassov et al. (2002). 



5. Monte Carlo Simulations 

As discussed in the previous section, the spectra of most sources are consistent with 
an absorbed power-law model. However, since the x^ statistic is not sensitive to small but 
systematic deviations that would result, for example, from the presence of a weak narrow line, 
the true null hypothesis probability may be lower than that inferred from the x^ statistic. 
In addition, the presence of a weak line may be overlooked when the spectra are binned to 
contain at least 20 counts per bin so that the x^ statistic can be used for parameter fitting. 

In this section, we describe the Monte Carlo (MC) simulations that are performed to 
test whether any of the discrete, unresolved line-like features are statistically significant. 
The method is identical to that adopted by Rutledge & Sako (2003), which was used to 
estimate the multi-trial statistical significance of the reported soft X-ray features observed 
in the spectrum of GRB011211 (Reeves et al. 2002). We briefly summarize the procedures 
below. 

For each source spectrum and time interval, we perform 10^ simulations based on the re- 
spective best-fit models. We denote the detected number of photons (source plus background 
events that lie within the extraction region and time selection) by Ctot and the estimated 
number of background photons by Cbkg, where Ctot > C'bkg- Since we do not know the true 
number of source photons and hence the source flux, we simulate a fixed number of source 
plus background photons. The number of background photons is estimated using a Poisson 
deviate with an average of Cbkg counts^^. Background events are randomly selected from 



^"^There are several alternative approaches one might adopt in choosing the number of counts to generate 
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analytic background spectral models, which are derived from fits to the spectra obtained 
from off- axis regions of the same data sets ( Chandra and XMM-Newton) or from long blank 
field observations {ASCA and Beppo-SAX). 

We generate counts spectra Ij for each simulation and convolve it through a matched 
filter to compute the C{Ei) values defined as, 



j[Ei-3a(Ei 



C{E, 
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where Ei is the energy in keV of PI channel i and cr{Ei) is the energy-dependent gaussian 
width of the instrument response kernel (see, Rutledge & Sako 2003). For each instrument, 
the shape of the response kernel is extracted at five energies - 0.5, 1, 2, 5, and 10 keV - from 
their respective response matrices and fit to a gaussian model. The best-fit cr{E) at these 
five points are then fit to a quadratic function of energy of the form. 



a{E) = Co + CiE, 



keV 



C2-EkeV> 



(3) 



where i?kcv is the energy in keV. The values for q derived in this manner for each of the 
instrument are listed in Table 4, and they reproduce the energy dependence of a to within 
~ 20% across the entire bandpass. The response kernels, in general, are asymmetric and 
non-gaussian. However, a gaussian reproduces the core of the kernel fairly well, and the 
results are not sensitive to their exact shapes outside of the core. For the HEG and MEG, 
we adopt a cubic polynomial of the form. 



O-(A) = Co + CiA^ + CsA^ + CgA^, 



(4) 



where the Ax is the wavelength in A. We adopt the parameters Ci as listed on the MIT 
HETG web site^^. We list them in Table 4 for completeness. Finally, for the LEG, we assume 
a constant wavelength dispersion of a = 0.0195 A, which is a good approximation in the 
2 — 20 A region. 



in each MC simulation. One possibility is to use a Poisson deviate with an average of Ctot and Cbkg counts 
for the total and background portion of the spectra, respectively. This approach, however, assumes that the 
observed flux is in fact the true source flux, which is, in general, not the case. We, therefore, choose to fix 
the total number of simulated counts to the observed value, since this is the only true observable parameter. 
Adopting a Poisson deviate of the observed number of counts would, in any case, result in a lower significance 
of any discrete feature, so our derived significances assuming a fixed number of counts can be regarded as 
upper limits. 
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http://space.mit.edu/HETG/technotes/ede_0102/ede_0102.html 
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As emphasized by Rutledge & Sako (2003), this procedure maximizes the signal-to- 
noise of hne-hke features with intrinsic widths that are smaller than or comparable to the 
instrument response. For spectra acquired with the CCD, the search is sensitive to lines 
with widths narrower than v ~ 10,000 km s~^. For the handful of high-resolution spectra 
in our sample, however, the velocity range spanned by adopting the instrument width is 
rather limited {v ~ 1000 km s^ ). In these cases, we use widths that are up to 16 times the 
instrument width to search for potentially broader lines as well. 

For each spectral bin, we then sort the 10^ C{Ei) values in increasing order and draw the 
1000th, 100th, and 10th highest values as the 99%, 99.9%, and 99.99% single-trial significance 
levels, respectively, as a function of energy. Similarly, we take the 1000th, 100th, and 10th 
lowest points, which represent 99%, 99.9%, and 99.99% lower limits, to identify possible 
absorption features in the spectra^^. The data are convolved in the exact same manner as 
the simulated events, and we compare the values of each energy bin to the upper and lower 
bounds. Figures 33 - 84 show the results for each source. 

We search through the matched-filtered spectrum for line-like features that exceed or fall 
below a certain confidence level. If a significant feature is found, the single-trial significance 
at the peak of the feature is recorded. We then return to the MC simulations and count the 
ones that contain one or more features that meet or exceed that same single-trial limit at 
any arbitrary energy. The number of such simulations are recorded and divided by the total 
number of simulations performed to yield the true multi-trial significance. Qualitatively, 
the multi-trial probability of detecting a line is higher by a factor of approximately the 
total number of resolution elements in the detector bandpass. For example, an emission 
feature with a single-trial significance of 99.99% (i.e., occurs once in 10'^ random trials at the 
specified energy) in a detector with 100 resolution elements is seen approximately once in 
~ 100 random trials at an arbitrary energy. The energy dependence of the resolving power 
and the effective area, however, complicates this and the simplest way to address this is 
through MC simulations. 



^^Note that, particularly for the 99.99% confidence limits, the MC function C{Ei) can be Poisson limited 
(that is, the few number of simulations produce strong point-to-point variations). When the data approach 
the confidence limit where the these variations can be seen in the simulations, we increased the number of 
simulations. In selected parts of the spectrum where the the number of counts per PI bin is extremely low, 
the C{Ei) can only possess a small set of discrete values, which results in wiggles in the confidence curves 
(see, e.g., the 99.99% curves in the PN spectra above ~ 2 keV, as well as those of the Chandra grating 
spectra) that cannot be smoothed out by increasing the number of simulations. The data, however, are 
convolved with the matched-filter in the exact same manner as described above, so the derived significances 
are robust. 
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This method provides a robust estimate of the statistical significances of discrete, unre- 
solved features observed in the spectra, assuming the energy of the feature is not known a 
priori. Reeves et al. (2003), Watson et al. (2002a) and Watson et al. (2003a) adopt a different 
method, where each of the simulated spectra are fit by models with and without emission 
lines. The changes in x^ are, then, recorded and compared to the observed Ax^ derived from 
the data. This method should, in principle, yield the probability that the observed features 
are due to chance coincidences. In practice, however, the significances derived using this 
method will be overestimated since, as described by Rutledge & Sako (2003), an automated 
X^-minimization is almost never capable of finding a true global minimum when the x^ sur- 
face is complicated and contain multiple local minima, for example, when the model involves 
multiple emission lines. Therefore, the inferred Ax^ between the models with and without 
emission lines for each simulation is underestimated. This is why the statistical significance 
derived by Reeves et al. (2003) for the features in the GRB011211 spectrum is significantly 
higher compared to that of Rutledge & Sako (2003). 

Finally, we make explicit our conventions for quoting an "equivalent gaussian sigma" 
for a given significance. When the MC simulation of M trials produces A^ values in excess of 
that observed at any given energy (single-trial), we give the probability of having observed 
such a high value produced from noise alone as Pexcess = N/M. This probability includes all 
phase-space for the observable, both above and below the median. Note that if our observed 
value is exactly at the median, then Pcxccss = 0.5. 

We then quote a "equivalent gaussian sigma" for this probability. This provides a short- 
hand and intuitive denotation of the significance of the observation. We adopt the convention 
that, for a probability Pexcess of producing a value at or above the observed value from noise 
alone, the equivalent gaussian sigma n is found from 

1 r°° 

Pexcess /t; — / 6 CtX, yO ) 
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which means that, for a measured value at the median of the simulated distribution (pexeess=0.5), 
this would be Oa in excess of the median value. We would call this feature to be signficant 
at the Oa level - and it would not be regarded as a detection of anything. In the case that 
Pexeess = 2 X 10~^, then n = — 4(T - or, a 4a deficit in counts beneath the median, which 
would be a significant detection of an absorption line-like feature. 

When multi-trial signficances as estimated by counting the N' simulations that show 
emission features with single-trial significances above a certain value, the equivalent gaussian 
sigma n' is derived from 

Pexcess ='W-/ C'^^'^^dx. (6) 
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which is a a one-sided distribution - i.e., if all of the MC simulations show an excess, then 
the significance is Oa. In comparison with Equ. 5, Equ. 6 gives a slightly higher value 
for the equilvalent gaussian sigma for the same Pexcess- Antonelli et al. (2000), Watson 
et al. (2002a), Butler et al. (2003b), and Watson et al. (2003a) have adopted Equ. 6 in 
deriving the equivalent gaussian sigma. The differences are relatively small especially at 
higher significances. 



6. Results of the Line Searches 

As shown in Figures 33 - 84, there are a number of discrete line-like features that exceed 
the 99.9% hmit, which approximately corresponds to 3. la at that particular energy (single 
trial) assuming a gaussian distribution of counts. In 26 cases listed in Table 5, the single-trial 
confidence exceeds 99.9% (3.1(t)^^. 

As we have noted earlier, we do not know a priori the expected energies of the lines, even 
in cases where the source redshift is known, since in some models the line emitting material 
can have an arbitrary velocity. In such cases, the relevant quantity is the chance probability 
of detecting any feature that exceeds the 99.9% limit (3.3a, multiple trials). As shown in 
the final column of Table 5, we find signals which meet this confidence limit in four datasets: 
in GRB011211 (at time t >5 ksec), GRB030227 (the complete dataset), GRB021004, and in 
GRB040106. 

Below, we first discuss our results in comparison with previous reports of line detection, 
followed by discussion of the four datasets with features with multi-trial probabilities with 
> 3.3(7 gaussian equivalent significance. 



6.1. Comparison with Previous Reports of Line Detection 

Below we compare our results to cases where line detections have been claimed in the 
published literature. 

GRB970508: Piro et al. (1999a) reported an unresolved line feature at 3.5 keV with a 
significance of 99.3% {2.7a) during time interval la during part of the afterglow. We find no 
evidence for a significant excess at 3.5 keV, or any energy (see Figure 33) during the same 



^^Note our use of nomenclature. We use the term "XX% confidence" , when XX% of the signals produced 
by background alone are at or below that detected. 
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time interval. The single-trial confidence at 3.5 keV is ~ 60% (~ 0.25a). The corresponding 
multi-trial confidence limit is only ~ 15% (~ 0.2o"). 

GRB970828: We find no evidence of any single-energy feature similar to the one 
claimed by Yoshida et al. (1999) during time interval B. However, the reported feature was 
claimed to be resolved with a centroid oi E = b.OA'^Q'lf keV and a width of o" = 0.3llo'3| keV 
(Yoshida et al. 1999, 2001). We see an excess at a similar energy in the both the SIS and GIS 
data (Figure 38 and 37), with single-trial significances of 99.2% (2.4(t) and 99.4% {2.5a), 
respectively. Our estimate of the multi-trial significance is only ~ 70% (~ l.Oo") in each 
instrument. Searches for lines with larger widths resulted in those with lower significance. It 
is interesting, however, that the residuals of the power-law model fit visually show excesses 
in all four detectors (see Figure 4). We also note that the S/N in this spectral range is only of 
order unity, which implies that the estimated significance is highly-dependent on the exact 
shape and magnitude of the local background spectrum. The reality of this excess feature 
should, therefore, be taken with caution. 

GRB991216: In the case of the GRB991216 Chandra spectrum, where a 4.7a detection 
of an emission line at -E ~ 3.5 keV has been claimed (Piro et al. 2000), we instead find only a 
99.0% (~ 2.3a) confidence single-trial fiuctuation in the dispersed HEG spectrum at a width 
of 16 times it's instrument resolution or FWHM ~ 0.7 keV (see Figure 40), which is close to 
value quoted in Piro et al. (2000). This fiuctuation is at a slightly higher energy (~ 3.8 keV). 
Our multi-trial confidence limit of this feature is 20% (0.3o"). We do not find any feature 
near 4.4 ± 0.5 keV, as reported as a marginal detection by Piro et al. (2000). The zeroth 
order spectrum appears to show an excess feature centered at -E ~ 3.4 keV, which is below 
that seen in the dispersed spectrum, (see 39), but is not significant. Searches for significant 
features that are narrower than 16 times the HEG resolution resulted in a non-detection. 

At least part of the source for the discrepancy between the results of our analysis and 
that of Piro et al. (2000) is the adopted continuum level. Piro et al. (2000) use a power- 
law slope and fiux determined from the zeroth order spectrum. However, as can be seen 
in their Figure 1, this continuum systematically underestimates the fiux across the entire 
HETG bandpass. Any fluctuation is, therefore, in excess of the adopted continuum flux. 
If we allow for a ~ 20% calibration uncertainty between the zeroth and flrst orders, the 
spectra of both orders are well-represented by a single continuum and, as a consequence, 
the single-trial signiflcance of the line decreases to 2.3a at a single energy. Inspection of the 
GRB020813 and GRB021004 spectra, both of which were collected in the same instrument 
conflguration, also show a consistent systematic ~ 20% discrepancy in the normalizations 
between the zeroth and flrst orders. 

Incidentally, we note that a feature seen in the MEG at E ~ 1.32 keV is of relatively 
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higher significance than the E ~ 3.5 keV feature reported by Piro et al. (2000), ahhough it 
is still extremely marginal. The formal significance is 99.91% single trial (~ 3.1a) and 89.2% 
multi-trial confidence (~ I.60"). 

GRB000214: Antonelli et al. (2000) find a feature at E ~ 4.7 keV feature which they 
claim, using an F-test to be significant at 99.73% at a single energy. In our simulations, 
we find a broad excess at similar energy (Figure 43), significant at ~ 99.92% (~ 3.2o") at a 
single energy. This is slightly higher than found by Antonelli et al. (2000). However, since 
the redshift is not known, the appropriate multi-trial significance is 99.03% (~ 2. 60"). 

GRB001025: Watson et al. (2002a) find a plasma model to be a better fit to the 
spectrum of this GRB than an absorbed power law, due to a soft excess observed between 
~ 0.5 and 2 keV. We find no significant excess at these energies (Figures 10, 45, and 46). 
The single-trial probabilities are all < 99.9% confidence, which is below that found with 
identical instrumentation for GRBO 11211 (see below). The multi-trial significance of the 
most significant feature (~ 0.7 keV excess seen in the PN) is ~ 65%. 

GRB010220: Watson et al. (2002a) claim the existence of a feature in the PN spectrum 
of this afterglow at 3.9 keV. As shown in Figures 11 and 47, we find no evidence for features 
at any energy being detected with single-energy significance >99% confidence {2.3a). The 
multi-trial significance limit would be, consequently, less than that found for GRB001025 
(see above), which is < 65% confidence (< 0.9o"). 

GRB011211: In GRB011211 the presence of several soft X-ray emission lines from 
mid-Z elements during the first 5 ksec of the observation has been claimed (Reeves et al. 
2003, 2002). We find one feature in the PN spectrum at E = 0.85 keV with 99.9% single-trial 
significance (see Figure 49), coincident with a line Reeves et al. (2002) identified with 97% 
confidence. This is consistent with the results of Rutledge & Sako (2003) and Borozdin & 
Trudolyubov (2003). The chance probability of finding such a feature at an arbitrary energy 
is ~ 15% (l.Aa). 

In the spectra acquired during the latter (t > 5 ksec) time interval, there is an excess 
feature in the PN spectrum (Figure 51) near the instrumental oxygen edge just above E ~ 
0.5 keV. This feature is statistically significant at the > 99.99% level at an arbitrary energy, 
so the formal confidence is > 3.89a. The feature is not seen in the MOS spectrum, which as 
an exposure (in cm^ s) a factor ~ 3 smaller at this energy. We cannot exclude the possibility 
that the excess is due to calibration uncertainty near the oxygen edge, although no such 
features have been reported in other high SNR data, nor did we find an excess at this energy 
in our analysis of the continuum source 3C273 (see below). What makes this particular case 
somewhat suspicious is that the background spectrum shows an excess similar to the one 
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seen in the GRB spectrum. 

GRB020813: Butler et al. (2003b) have reported evidence for a hue at 1.3 keV. We do 
find a narrow fiuctuation at approximately this energy in the HEG spectrum at ~ 99.99% 
confidence (Figure 60), single trial. The multi-trial confidence is, however, only 97.86% (2.3 
a). Note that the decrease in confidences between single-trial and multi-trial is greater for 
grating observations, due to the larger number of independent resolution elements in grating 
spectra. 

Inspection of the raw counts in each of the dispersion orders shows an excess in only the 
m = +1 order of the HEG. Of the four spectral orders that are considered (m = ±1 in HEG 
and MEG), this is the one with the lowest effective area in this spectral range. The feature is 
not detected in any of the other spectral orders including the zeroth order. In particular, the 
effective area of the MEG m = —1 order is ~ 4 times larger than that of the HEG m = +1, 
which makes the reality of this detection somewhat suspicious. 

GRB030227: This is the second highest signal-to-noise observation in our sample only 
after GRB040106. Watson et al. (2003a) find four emission lines between 0.62 and 1.67 keV 
in the last 11 ksec of the PN spectrum of this event. In our analysis (Figures 66 - 72), we 
find no significant (at > 99.9%) features in the spectrum during this time interval (a feature 
near 0.84 keV is close to 99.9% confidence, single trial). The most significant feature, with 
single-trial detection confidence of < 99.9%, has multi-trial detection confidence of < 85% 
(< lAa). 

We do, however, find several fiuctuations that are in excess of the 99.99% limits in the 
cumulative PN and MOS spectra (Figures 66 and 69). One of these is in the 1.60 — 1.69 keV 
range, coincident with a feature identified previously by Mereghetti et al. (2003). However, 
the multi-trial significance is marginal (98.93% confidence detection; 2.60"). Two features at 
3.95 — 4.06 keV and 4.33 — 4.39 keV have multi-trial significance greater than 3.3a (3.54o" 
and 3.62o", respectively). The former is particularly interesting, since there is a feature at 
approximately the same energy {E ~ 3.9 keV) in the MOS as well, just above the 99.9% 
limit. 

The background spectra of the PN or the MOS both do not show discrete features near 
this energy. To check for other possible calibration uncertainties, we have analyzed the 
XMM-Newton data of a featureless continuum source, 3C273, to see whether excesses of 
similar significance do indeed appear in the spectra of other sources. We have extracted only 
a small fraction of the events so that the spectrum contains the same number of counts in the 
1-5 keV range as in the GRB030227 spectrum. We fit the data to an absorbed power-law 
model and performed MC simulations to derive formal significances of the local fiuctuations. 
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Interestingly, no features in excess of 99.99% single-trial confidence were detected in either 
the PN or the MOS data of 3C273. 

In summary, in most cases where lines have been claimed, we either find the fluctuations 
to be of lower significance, or to be absent. A 3ct feature at a given energy in a spectrum 
with few counts visually appears more significant (i.e., higher equivalent width) than in a 
spectrum of high statistical quality. These fluctuations are also transient and would appear 
only during a segment of the total exposure. This trend is similar to the ones seen in 
many observations in the sense that (1) the "lines" are transient, (2) almost always ~ 3a 
irrespective of the statistical quality of the data, and (3) high equivalent width. 



6.2. Other Sources that Exhibit Line-like Fluctuations 

As summarized in Table 5, several sources other than those reported in the literature 
show features in excess/deficit of the 99.9% single-trial upper/lower limits. Six of these 
features in four spectra meet a 3a multi-trial signiflcance cut-off: 

1. GRB011211: As discussed in §6.1, this source shows a highly signiflcant excess in 
the 0.51 — 0.64 keV energy range, with a multi-trial signiflcance in excess of 3.89a. If 
modeled as a gaussian line, the feature is unresolved with a centroid of £^ = 0.58 keV. 
The observed line flux and equivalent width [EW) in the PN data is F = 1.1 x 
10^^ photons cm~^ s"^ and EW = 70 eV, respectively. This is one of the two most 
statistically signiflcant excess in the present analysis. The la upper limit in the MOS 
data is EW = 27 eV. The rest frame energy for this feature is -E = 1.6 — 2.0 keV, which 
corresponds to a K-shell transition array of neutral to hydrogen-like Si. The excess 
at this same energy in the background spectrum (discussed above), however, makes 
this case somewhat suspicious. We leave further analysis with a reflned calibration to 
future work. 

2. GRB021004: One feature in the zeroth order spectrum at E = 1.60 — 1.69 keV is 
signiflcant at > 99.99% (> 3.72a) single-trial and at 99.92% (3.35a) multi-trial. This 
feature lies below the lower 99.99% limit, which implies an absorption feature. The 
line is again unresolved, with a centroid oi E = 1.64 keV. The absorption line flux 
and equivalent widths are F = —1.0 x 10^^ photons cm^^ s~^ and EW = —91 eV, 
respectively. The rest frame energy oi E = 5.3 — 5.6 keV does not correspond to any 
strong atomic transition. An obvious feature at the same energy is not observed in 
the dispersed high resolution spectrum, with a la lower limit of —26 eV, but is still 
consistent with the zeroth order spectrum. 
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3. GRB030227: We find six features in tlie PN spectrum and two in tlie MOS spectrum 
wliicli liave single trial confidence levels >99.99%; two of these have multi-trial confi- 
dences levels > 3o" observed in the PN detector. The centroids are at E = 4.00 keV 
and E = 4.36 keV, with emission line fluxes oi F = 1.0 x 10"^ photons cm~^ s~^ and 
F = 1.1 X 10~^ photons cm~^ s~^, respectively. The corresponding equivalent widths 
are EW = 66 eV and EW = 83 eV. The upper limits derived from the MOS spectra 
are EW = 93 eV and EW = 32 eV for the 4.00 keV and 4.36 keV features, respectively. 
The host redshift is not known and so the features cannot be identified with known 
atomic transitions. 

4. GRB040106: Three features in the PN and one in the MOS have single-trial signifi- 
cances in excess of 99.9%. The feature at ~ 0.6 keV is formally significant at > 3.89 a, 
and its structure is similar to the one detected in the t > 5 ksec of GRB011211. The 
corresponding flux observed in the PN detector is F = 1.4 x 10~^ photons cm"^ s~^ 
{EW = 39 eV) with a line centroid oi E = 0.66 keV. The upper limit derived from the 
MOS data is F = 9.1 x 10"^ photons cm^^ g-^ {EW = 26 eV). The ~ 1.33 keV feature 
is an absorption line-like feature with a multi-trial confldence of 99.89% or 3.26o". The 
absorption flux observed in the PN spectrum is F = —2.2 x 10~^ photons cm~^ s~^ 
{EW = —18 eV), and the lower limit in the MOS data is F = 9x 10"'' photons cm"^ s~^ 
{EW = 7 eV). Once gain, since no redshift measurement of the host is available to 
date, we cannot identify any of the features with strong atomic transitions. 



6.3. Measurements of Equivalent Hydrogen Column Density in the Host 

Galama & Wijers (2001) analyzed eight GRB X-ray afterglows, of which only one of 
them (GRB970508) is part of our sample, and found that they show evidence for absorption 
in excess of the Galactic line of sight values. If the absorbing material is at the redshift of 
the GRB, the columns are large, so that combined with low Ay as measured from multi-color 
optical data, Galama & Wijers (2001) take this to be evidence for dense ionized material 
surrounding the progenitor. The low Ay values can be accounted for by assuming the 
material is ionized by radiation from the explosion. 

If we consider all of our measurements, 14 out of the 21 total have a host column with 
either a best-flt value of zero, or within 2a of zero (Table 3). In four cases, GRB970508, 
GRB020322, GRB020405, and GRB030227, there is strong evidence for non-zero host Nh 
at > 4o", where the error includes only statistical uncertainty, but does not include any 
estimate of the error in the assumed Galactic Nh- Of the three most signiflcant (GRB020322, 
GRB030227, and GRB031203) the host redshift of two (GRB020322 and GRB030227) is. 
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unfortunately, unknown (we assume z = 1), so the measurement is problematic. GRB031203 
does appear to exhibit a column density that is slightly higher, though statistically significant, 
than that of Dickey & Lockman (1990), but it is consistent with the dust echo measurements 
as described in Vaughan et al. (2003). Therefore, in a majority of cases, we find no evidence 
for significant absorbing columns local to the burst. Of course, columns less than a few times 
10^^ cm~^ are difficult to measure given the significant redshifts for the events. 

Stratta et al. (2004) examined 13 afterglows observed with Beppo-SAX, which includes 
the two sources presented in this paper, and found that 2 sources show highly-significant 
absorption above the Galactic value. Both GRB970508 and GRB000214 do not appear to 
show any obvious signatures of substantial host galaxy absorption, which is consistent with 
what we find. The upper-limits they derive for the remaining 11 sources are comparable to 
or above the magnitude of the detections, as are the upper-limits we derive, implying that 
the non-detections are due to instrumental sensitivity. 



6.4. Consistency with Prior Detection Claims 

Our analysis at first glance appears to be inconsistent with previous work, in that we 
do not confirm any prior claimed detection to be significant. In addition, the most signif- 
icant deviations we find have not been previously reported, even in spectra with published 
analyses. In some cases we find deviations at the reported energy, but at statistical signif- 
icance substantially lower than claimed. In a few cases we do not see positive fiuctuations 
at the claimed energy, and in three cases we find fiuctuations at > 3o" that have not been 
reported in spectra analyzed by others. An important point to note is that our analyses, by 
no means, rule out the possibility of discrete features in any of the spectra. The presence of 
weak features is certainly possible, but, in most cases, they are not required by the data in 
a model-independent way. 

In most cases we can understand the origin of the discrepancy. In GRB991216 we 
disagree with the continuum level adopted by Piro et al. (2000). If the continuum model 
lies below most of the data points (see Figure 1 of Piro et al. 2000), the inferred significance 
of the line would consequently be overestimated, although we find that renormalizing the 
continuum to the one adopted by Piro et al. (2000) yields a single-trial significance of ~ 3.2o", 
which is still lower than claimed by those authors (~ 4.7o"). In most other cases, we have 
probably adopted different background regions for subtraction, so that if reported features 
are background fluctuations, we would expect to find them at different energy, or at different 
level depending on the exact background region used. In any case, a robust detection should 
not be highly-sensitive to the exact nature of the assumed background spectrum. 
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Our analysis results are in rough agreement with most of those published recently by 
Butler et al. (2004), who have also re-analyzed a small subset of the data presented here using 
a similar statistical approach. There are, however, notable differences in the quoted multi- 
trial significances of the features in the high-resolution grating spectra of GRB991216 and 
GRB020813. They find that the claimed feature in the first order spectrum of GRB991216 
(Piro et al. 2000) is significant at > 3.7cr multi-trial, while our analysis yields a significance 
of at most ~ I.Sct (see, also. Figure 40 and 41). Similarly, they find that the significance of 
the line claimed in GRB020813 by Butler et al. (2003b) is 3.5cr multi-trial, while we find it 
to be significant at only 2.1a. These discrepancies are most likely due to the fact that Butler 
et al. (2004) choose to simulate counts based on their spectral fit to unbinned spectra, which 
predicts a count rate that is lower than the observed value. This means that fluctucation 
present in the spectrum appear preferentially as positive emission line-like fluctuations. We 
believe that our approach of fixing the number of simulated counts to the obseved value is 
more robust, unless there are good reasons to believe a priori that the true continuum level 
is lower than what one infers from a blind fit. 

We point out that a true multi-trial significance should also account for the total number 
of spectra inspected, and depends on whether the observation was divided into several time 
segments or if data from more than one detector with comparable S/N and spectral coverage 
were available. It appears that none of the authors have properly derived the significances 
of the features that are seen only in selected time segments of a longer observation and in 
a single detector. For example, if a 3(T line is detected only in the PN spectrum aqcuired 
in the first 5 ksec of a 40 ksec exposure, then the true multi-trial probability that this is 
not due to a statistical fluctuation is ~ 2 x 40/8 = 10 times lower or ~ 2.2o" (the factor of 
~ 2 comes from the fact that the two MOS cameras combined provide data of comparable 
sensitivity). If the S/N for each time segment and/or detector are not identical, one must 
carefully weigh each trial by estimating its sensitivity to the detection of the line seen in the 
other dataset(s). 

In all cases, the authors have not explicitly stated all of the time segments used in their 
line searches. As stated earlier, none of the features are observed in more than one detector 
and this has not been factored into their derived significances. We have not factored these 
criteria into the quoted significance for the features that we detect either, since we make 
no claims as to the reality, interpretation, or ubiquity of features such as the ones that 
appear in GRB011211, GRB021004, GRB030227, or GRB040106. Rather, we find these to 
be motivation for significantly higher signal-to-noise observations, and we call into question 
prior statements regarding the detection of features with a specific interpretation. 
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7. Summary and Conclusions 

We have performed detailed, uniform analyses of 21 of the brightest GRB X-ray af- 
terglows, and have derived the statistical significances of line-like fluctuations that appear 
in the spectra. We find that in most cases, the data are well-fit by an absorbed power- 
law spectral model, with the significance of any deviations being marginal. This includes 
several cases where either lines, or deviations from a power-law continuum have previ- 
ously been reported at reasonable significance. In particular, this includes GRB970508, 
GRB970828, GRB991216, and GRB000214 where features from Fe emission have been re- 
ported, GRB011211, GRB010220, and GRB030227 where a thermal emission hues from low-Z 
elements have been claimed, and GRB001025, where one prior analysis found the spectrum 
to be better described by a thermal emission model than by a power law. 

There are four interesting exceptions that we find are worth attention; the cumulative 
spectrum of GRB030227, the t > 5 ksec spectrum of GRB011211, the zeroth order spectrum of 
GRB021004, and GRB040106. In addition, there is one interesting case (GRB970828), where 
a broad excess is visible in all four detectors on ASCA, although the formal significance is 
marginal and depends highly on the assumed background spectrum. 

The absorption line-like feature seen in the zeroth order spectrum of GRB021004 is 
formally significant at 3.35o" multi-trial. This feature corresponds to an energy of ~ 5.5 keV 
in the rest frame of the GRB, and does not have an obvious identification. It is not observed 
in the dispersed spectrum. We also note that there is a local minimum in the effective area 
at this energy, which makes the feature appear more significant visually, as in the case for 
GRB991216, where the claimed emission line lies at a local maximum of the efficiency curve. 

The t > 5 ksec interval of the GRB011211 spectrum and the cumulative spectrum of 
GRB040106 both exhibit features at -E ~ 0.5 keV. These are by far the most statistically 
significant in our sample with multi-trial significances exceeding 3.89(7. The host's known 
redshift of GRB011211 places the feature at a rest frame of i? ~ 1.7 keV, which may be 
identified as Si K fiuorescence, but could easily be due to other elements moving at arbitrary 
velocities. As noted above, the observed energies correspond also to a local fluctuation in 
the background spectrum as well as an instrumental absorption edge of atomic or molecular 
oxygen. 

In GRB030227, Watson et al. (2003a) report the detection of multiple emission lines 
during the last ~ 1 1 ksec of the observation. We do not conflrm this result and we flnd that 
the strongest feature at -E ~ 0.85 keV is signiflcant at < 99.9% for a single trial. However, 
the cumulative spectrum of the entire observation does appear to show fluctuations with 
substantially higher signiflcance than those reported by Watson et al. (2003a). One of these. 
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at i? ~ 4 keV, is significant at 3.540". Unfortunately tliere is no measured redsliift for tliis 
event, so tlie rest-frame energy is unknown. 

Although our analysis rules out all of the reported features, we do find some evidence, 
in four cases, for marginally significant spectral features at energies that do not suggest a 
clear interpretation. We cannot rule out the possibility that relatively weak line features 
are present in the X-ray spectra of some, if no all, GRBs. Given the diagnostic power of 
such detections, longer exposures of bright events are certainly warranted. Our work does, 
however, call into question the reality of all of the features detected to-date, as well as the 
interpretation of these features as either iron, or specific low-Z elements. 

Swift, which will observe a large number of GRBs early on, when they are bright, has 
the potential and opportunity to confirm the existence of X-ray line emission. Although the 
effective area of the Swift X-ray Telescope (XRT^^) is a factor of ~ 10 lower than that of the 
XMM-Newton PN camera and a factor of ~ 20 lower than the combined EPIC effective area, 
it is capable of collecting data starting only minutes after the burst compared to ~ 0.5 days 
for both XMM-Newton and Chandra. For a flux decay time slope of a = 1 — 2, where the 
X-ray flux Fx oc t^", a typical burst will be 2 — 3 orders of magnitude brighter during the 
beginning of a Swift observation and the XRT data on average will be of higher statistical 
quality by roughly an order of magnitude. The late afterglow > 0.5 days after the burst will 
still be better studied with the CCDs on XMM-Newton, Chandra, as well as ASTR0-E2. If, 
however, the lines are narrow [v < 1000 km s"^), the Chandra gratings, the XMM-Newton 
RGS, and the X-ray Spectrometer on ASTR0-E2 wiW provide more sensitive measurements on 
most bursts. Clearly, in any case, this is not a straightforward measurement, and a carefully 
planned observation strategy is required to maximize the detection probability. Monitoring 
of the early X-ray flux with Swift and follow-up observations of the brightest afterglows with 
XMM-Newton, Chandra, and ASTR0-E2m][ minimize the risk of obtaining an underexposed 
spectrum, although this does not by any means guarantee the detection of lines. Of course, 
the ultimate goal is to take one step further and actually use X-ray spectroscopic diagnostics 
to infer important physical information about the burst. Note, however, that the physical 
conditions of the X-ray emission regions during the early afterglow (accessible only with 
Swift) can be vastly different from those approximately a day after the burst, and it is not 
clear whether a direct comparison of measurements will be meaningful. 

We finally point out that spectral identifications is also not a straightforward task. This 
is obvious when only a single significant feature is detected, but is also true even when 
multiple lines are detected in the spectrum. The energies of the hydrogen-like Lya lines 
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of neighboring abundant metals, for example, are spaced roughly equally and they can be 
misinterpreted or their identifications can at best be ambiguous. If a significant amount 
of iron is present, the spectrum may exhibit broad bumps that may dominate over narrow 
isolated lines from other elements and may not appear as prominent line emission. Although 
arbitrary bulk velocities can be invoked to match transitions of specific elements, any con- 
straint on the burst redshift (e.g., from optical absorption/emission line measurements) is 
extremely valuable. 

If the combination of Swift, XMM-Newton, Chandra, and ASTR0-E2 is unable to find 
similar features in spectra with large numbers of counts, then the ~ 3cr detections found here 
are likely spurious. If this is the case, tight upper limits on the line fiuxes and equivalent 
widths will provide invaluable constraints on the geometry and chemical abundances of the 
burst environment, which can then be used as quantitative tests of GRB afterglow models. 
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Table 1. Observations of Bright GRB X-ray Afterglows 



Burst 



Mission 



Instrument (s) 



T,t 



Exposure 



970508 


Beppo-SAX 


LEGS & MEGS 


0.4340 


15.4 


0.835 


970828 


ASCA 


SIS & GIS 


1.1906 


33.8 


0.9578 


991216 


Chandra 


HETG/AGIS-S 


1.5215 


9.65 


1.02 


000210 


Chandra 


AGIS-S 


0.7778 


8.93 


0.8463 


000214 


Beppo-SAX 


LEGS & MEGS 


0.5237 


43.5 




000926 


Chandra 


AGIS-S 


2.6861 


10.2 


2.066 


001025 


XMM-Newton 


EPIG-PN 


1.9092 


23.7 








EPIG-MOS 


1.8814 


37.1 




010220 


XMM-Newton 


EPIG-PN 


0.6247 


12.1 




011030 


Chandra 


AGIS-S 


10.4592 


46.6 




011211 


XMM-Newton 


EPIG-PN 


0.5029 


26.8 


2.140 






EPIG-MOS 


0.4851 


31.2 




020321 


XMM-Newton 


EPIG-PN 


0.4289 


36.3 








EPIG-MOS 


0.4119 


40.0 




020322 


XMM-Newton 


EPIG-PN 


0.6428 


23.8 








EPIG-MOS 


0.6222 


28.2 




020405 


Chandra 


LETG/AGIS-S 


1.6825 


50.6 


0.690 


020813 


Chandra 


HETG/AGIS-S 


0.7367 


76.9 


1.254 


021004 


Chandra 


HETG/AGIS-S 


0.8533 


86.7 


2.323 


030226 


Chandra 


AGIS-S 


1.5436 


36.4 


1.986 


030227 


XMM-Newton 


EPIG-PN 


0.3604 


33.3 








EPIG-MOS 


0.3440 


35.9 




030328 


Chandra 


LETG/AGIS-S 


0.6387 


92.7 


1.520 


030329 


XMM-Newton 


EPIG-PN 


37.0577 


29.6 


0.1685 






EPIG-MOS 


37.0423 


32.4 








EPIG-PN 


60.3962 


38.8 








EPIG-MOS 


60.3807 


46.5 




031203 


XMM-Newton 


EPIG-PN 


0.2723 


50.7 


0.105 






EPIG-MOS 


0.2567 


57.4 








EPIG-PN 


2.8549 


36.5 








EPIG-MOS 


2.8414 


34.4 




040106 


XMM-Newton 


EPIG-PN 


0.2552 


37.1 








EPIG-MOS 


0.2344 


42.6 





''start of the X-ray observation in days sinee the time of the GRB in the 
observers' frame 

''total exposure in ksec used in the analysis 

'^spectroscopic redshift of the host galaxy inferred from non-X-ray data 



Table 2. Summary of Reported Emission Lines in the Literature 



Burst Epoch/order Line Energy (kcV) Line Width (keV) Model/Identification Significance (cr)^ 



Reference 



970508 
970828 
991216 

000214 
001025 



la 

B 
\m\ = 1 
m = 



010220 
011211 (t < 



kscc) 



020813 \m\ = 1 

030227 (last 11 ksec) 



3.5 

3.49 ± 0.06 
4.4 ±0.5 
4.7 ±0.2 

1.16 ±0.05 
1.64 ±0.07 

2.2 ±0.1 

4 7+0.8 
^•'-0.4 

~3.9 
0.44 ± 0.04 
0.71 ±0.02 
0.88 ± 0.02 
1.22 ±0.03 
1.46 ±0.07 
1.31 ±0.01 

1-0.03 



0.62 



0.02 

+0.02 

0.03 



0.86: 
1.11 ±0.02 

1 ot;+0.04 
^■•-'°-0.03 

1.66 ±0.04 



1.67 



+0.01 
-0.03 



£0.5 


FeKa 


99.3% {2.7a) 


Piro ct al. {1999a) 


n O1+0.38 


Fe XXVI Lya or RRC*^ 


98.3% {2.4a) 


Yoshida et al. (1999, 2001) 


0.23 ±0.07 


Fe XXVI Lyo 


99.99974% (4.7a) 


Piro ct al. (2000) 




Fe XXVI RRC 


99.5% {2.8a) 


Piro ct al. (2000) 


unresolved 


Fe Kq 


99.85% (3.2a) 


Antonelli et al. (2000) 


0.10 ±0.03 


MgXII 


99.94% (3.4a) 


Watson et al. {2002a) 




Si XIV 


99.92% (3.4a) 






S XVI 


99.92% (3.4a) 






Ar XVIII 


98% (2.3a) 






Ni XXVIII 


88% (1.6a) 




unresolved 


Ni XXVIII Lyo 


99.0% (2.6a) 


Watson et al. (2002a) 


unresolved 


MgXI 
Si XIV 
S XVI 
Ar XVIII 
CaXX 


99.97% (3.6a)d 


Reeves et al. (2003) 


0.008 ± 0.004 


S XVI Lya 


99.90% (3.3a) 


Butler et al. (2003b) 


unresolved 


MgXII 


97% (2.2a) 


Watson et al. {2003a) 




Si XIV 


99.98% (3.8a) 






S XVI 


99.96% (3.5a) 






Ar XVII* or XVIII 


92% (1.7a) 






Ca XIX or XX 


99% (2.5a) 






Fe Ka 


99.86% (3.2a) 


Mereghetti et al. (2003) 



00 



''Probability of producing the observed data from noise only, as reported in the reference. The a value is the equivalent gaussian number of 
standard deviations corresponding to the stated confidence level. 

''Radiative recombination continuum 

"^The lines are collectively interpreted as due to thermal line emission at a rcdshift of 2 = 0.7_q'j. 

'^Combined multi-trial significance of the five lines, interpreted as due to thermal line emission at a rcdshift of z = 1.91 ± 0.06. 

"The five lines are interpreted as due to thermal line emission at a rcdshift of z = 1.39_Qgg. 

Correcting a typo in Watson et al. (2003a). 
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Table 3. Continuum Spectral Parameters 



Burst 


Epoch/order 


Ctot'^ 


Cfekg"^ 


^Galc 


^hostd 


pe 


i^lkcV 


f 




2 / 


P(X,I^)'' 


970508 


la 


137 


34 


0.051 


ioi« 


o 1+0.5 
•^■^-0.4 


9 fi+0.5 


X 


10" 


< 


2.6/5 


0.76 




lb 


188 


56 


0.051 


7tV 


2 n+10 


c:+0.4 
"^•^-0.3 


X 


10- 


1 


2.4/6 


0.88 


970828 


A+C 


1922 


953 


0.036 


i-5t;;? 


2 4+0-4 


2 0+1-0 


X 


10- 


1 


100/86 


0.14 




B 


1215 


515 


0.036 


o.o«:« 


r, 1+0.3 

^■^-0.3 


i-7t;];^6 


X 


10- 


1 


58/54 


0.33 


991216 


\m\ = 1 


762 


13 


0.20 


o.st^-j 


1 6+"-l 
-■-■0-0.1 


8.2tO-J 


X 


10- 


1 


153/130 






m = 


604 


4 


0.20 


o.o«:f, 


1 7+0.1 
-■-■'-0.1 


6 7+0-5 
"•'-0.5 


X 


10- 


1 


46/32 


0.05 


000210 






552 


2 


0.025 


3.7ti:i 


-I q+0-1 


8.810-? 


X 


10- 


'> 


26/21 


0.22 


000214 






378 


134 


0.058 


o.otj;^ 


9 0+0-3 
^■'J-0.3 


1 4-1-0.5 


X 


10- 


1 


28/18 


0.07 


000926 






267 


1 


0.027 


n o+o* 


-1 q+0-3 


9+0.4 


X 


10- 


'> 


15/13 


0.20 


001025 






1249 


313 


0.061 


n. iq+0.07 


2 fi+0-4 


9 c:+0.7 
^•"-0.5 


X 


10- 


J 


43/63 


0.98 


010220 






247 


124 


0.86 


o.itj;? 


1 R+1-0 
-■-■=-0.4 


1 9+0-^ 

-■-•^-0.6 


X 


10- 


'' 


6.5/9 


0.69 


011030 






376 


6 


0.10 


o.olU.i 


1 t;+0.3 
-■-■0-0.2 


l-^l 


X 


10- 


~i 


3.1/14 


1.00 


011211 


t < 5 kscc 


936 


58 


0.042 




nn + 0.2 
^■^-0.2 


7.8tJ;? 


X 


10- 


'' 


38/40 


0.56 




t > 5 kscc 


3046 


294 


0.042 


0.2tO-J 


9 o+O.l 
^■■^-0.2 


4 6+0-5 
^•"-0.1 


X 


10- 


J 


99/129 


0.98 


020321 




1664 


487 


0.082 


0.02t°«2 


n 1+0.2 
^■^-0.3 


9 9+2.3 


X 


10- 


s 


33/45 


0.91 


020322 




5791 


1252 


0.046 


20+°''3 


o O + O.l 
^■■'-0.1 


1 1+0-1 

1-1-0.1 


X 


10- 


1 


228/241 


0.72 


020405 


\m\ = 1 


1242 


84 


0.043 


o.4t»:l 


1 7+0.1 
-■-■'-0.1 


9 fl+0.1 
^•"-0.1 


X 


10- 


1 


1364/1318 






m = 


604 


4 


0.043 


o.6«:? 


2 0+0-2 

^■'J-o.i 


Y+0.8 
■^•'-0.7 


X 


10- 


1 


32/33 


0.49 


020813 


\m\ = 1 


4659 


117 


0.075 


o.olS;Si 


1 R+0-1 
-■-■=-0.1 


6 -^+0-2 
D.d_o.2 


X 


10- 


1 


455/348 






m = 


3679 


77 


0.075 


o.ot^.J 


1 s + "-i 


5 2+0-2 


X 


10- 


1 


145/142 


0.42 


021004 


\m\ = 1 


1634 


109 


0.011 


o.ot^.i 


O 1+0.1 

^■^-0.1 


9 f.+0.1 


X 


10- 


1 


440/348 






m = 


1252 


5 


0.011 


o.ot»:f, 


2 0+0-1 
^■'J-o.i 


1 6+01 

i-"-o.i 


X 


10- 


1 


77/52 


0.02 


030226 




345 


1 


0.016 


n 0+0-5 


9 1+0.2 
^■1-0.2 


i-3t'^;l 


X 


10- 


'' 


11/13 


0.59 


030227 


total 


13676 


2348 


0.22 


60+0'''* 


1 07+0.02 
^•^'-0.02 


2.34±0:r4 X 10- 


-4 


303/278 


0.12 




last 11 kscc 


2777 


171 


0.22 


0.7tS1 


-1 q+0-1 


1 0+0.2 

l-°-0.2 


X 


10- 


1 


115/124 


0.70 


030328 


\m\ = 1 


562 


131 


0.043 


o.ot'^;f, 


1 q+0-1 
-■-■^-0.1 


e.ei'clj 


X 


10- 


) 


1172/1318 






m = 


529 


18 


0.043 


o-ist'^i^ 


2 1+0-2 
^■1-0.1 


7.^t\l 


X 


10- 


' 


30/21 


0.09 


030329 


I 


665 


80 


0.021 


00+"''^ 
'^^ 'J'-' -0.00 


2 n+0-2 

^■^-0.2 


7 4+0.8 
'•^-0.9 


X 


10- 


3 


37/37 


0.48 




II 


543 


143 


0.021 


02+0''7 


2 0+0-3 
^■'J-0.2 


3.9+0-0 


X 


10- 


s 


34/31 


0.32 


031203 


I 


7012 


647 


0.59 


n r,q + 0.04 

'J-^y-0.03 


1 S+0-1 

^■«-o.i 


1-itri 


X 


10- 


1 


218/208 


0.32 




II 


1194 


381 


0.59 


0.80_0 ;^g 


1 q+0-1 
1-^-0. 1 


5-8t21 


X 


10- 


J 


41/50 


0.81 


040106 




23190 


1420 


0.086 


00+"''2 
'^^ 'J'-' -0.00 


1 c-f,+0.02 

i-ou_o.oi 


1.8810.0; 


X 10- 


-4 


622/603 


0.29 



^total number of source and background counts detected within the source extraction region 
''estimated number of background counts within the extraction region 
■^Galactic column density in multiples of 1022 cm~2 from Dickey & Lockman (1990) 
"^inferred column density at the redshift of the host galaxy in multiples of 1022 cm-2 
"photon index 



'continuum flux at 1 keV in units of ph cm 
SnuU hypothesis probability 



keV- 
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Table 4. Parameters for the Matched Fiher KerneP 



Instrument 


Co 


Cl 


C2 


C3 


SIS 


4.43 X 10^2 


1.28 X 10-2 


2.52 X 10-5 




GIS 


4.13 X 10^2 


3.82 X 10-2 


-1.60 X 10-3 




MEGS 


5.20 X 10-2 


3.32 X 10-2 


-1.12 X 10-3 




AGIS 


3.96 X 10-2 


5.28 X 10-3 


-1.31 X 10-^ 




MOS 


2.50 X 10-2 


8.77 X 10-3 


-2.79 X 10-^ 




PN 


3.42 X 10-2 


6.63 X 10-3 


-1.46 X 10-^ 




HEG 


1.00 X 10-2 


-1.41 X 10-^ 


1.46 X 10-5 


-1.963 X 10-^ 


MEG 


1.89 X 10-2 


-3.39 X 10-^ 


2.98 X 10-5 


-4.110 X 10-^ 


LEG 


1.95 X 10-2 









'^The Q represent coefficients for the i order term in the polynomial for 
(t{E) in keV and o"(A) in A for non-dispersive and dispersive instruments, 
respectively. See Equ. 3 and 4 for details. 



-46- 



Table 5. Summary of Single- and Multi-trial Significances of Features with Single- Trial 

Significances >99.9%'' 



Burst 


Instrument 


Epoch/order 


Energy Range (keV)'' 


Single-trial 


Multi-trial 


991216 


MEG 


mi = 1 


1.317- 


- 1.326 


99.98% 


89.20% 


000214 


MEGS 




4.60- 


-4.80 


99.92% 


99.03% 


001025 


MOS 




2.60- 


-2.63 


> 99.99% 


99.23% 


011211 


PN 


t < 5 ksec 


0.88- 


-0.91 


99.96% 


89.09% 




PN 


t > 5 ksec 


0.51 - 


- 0.64 


> 99.99% 


> 99.99% 


020321 


MOS 




1.01 - 


- 1.03 


99.95% 


94.85% 




MOS 




7.58- 


-7.66 


> 99.99% 


99.36% 


020322 


PN 




6.15- 


-6.18 


99.93% 


89.15% 




MOS 




2.09- 


- 2.12 


99.92% 


80.13% 




MOS 




4.97- 


-5.08 


99.90% 


72.26% 


020405 


LETG 


mi = 1 


0.733 - 


- 0.740 


99.97% 


98.23% 




LETG 


m = 1 


0.941 - 


- 0.953 


99.96% 


96.98% 


020813 


HEG 


m = 1 


1.303 - 


- 1.322 


> 99.99% 


99.40% 


021004 


ACIS-S 


mi = 


1.60- 


- 1.69 


> 99.99% 


99.92% 


030227 


PN 


total 


1.51 - 


- 1.55 


> 99.99% 


98.92% 




PN 


total 


1.66- 


- 1.69 


> 99.99% 


98.93% 




PN 


total 


1.80- 


- 1.96 


> 99.99% 


98.99% 




PN 


total 


2.03- 


- 2.11 


> 99.99% 


98.93% 




PN 


total 


3.95- 


-4.06 


> 99.99% 


99.96% 




PN 


total 


4.33- 


-4.39 


> 99.99% 


99.97% 




MOS 


total 


1.85- 


- 1.88 


> 99.99% 


99.25% 




MOS 


total 


2.62 - 


-2.66 


> 99.99% 


99.50% 


040106 


PN 




0.54- 


-0.70 


> 99.99% 


> 99.99% 




PN 




1.30- 


- 1.37 


> 99.99% 


99.89% 




PN 




1.95- 


- 1.96 


> 99.99% 


98.58% 




PN 




4.07- 


-4.11 


99.94% 


93.49% 



''Only discrete features that are in excess of 99.9% single-trial significance (3.29o") are listed 
with the exception of GRB030227, where we adopt a lower cut-off of 99.99% {3.89(t) owing to 
the relatively large number of features that lie beyond the 99.9% upper and lower limits. 

''For features with single trial significances between 99.9% and 99.99%, this corresponds to the 
range in energy where the data exceed the 99.9% limits. For features with single trial significances 
above 99.99%, this corresponds to the range in energy where the data exceed the 99.99% limits 
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1 1 r 



1 ' I I 

GRB970508 ■ 
Beppo-SAX (la) 



A LEGS 
• MEGS 




1 2 5 

channel energy (keV) 

Fig. 1. — The observed X-ray spectrum of GRB970508 time segment la. In this figure and the 
similar ones that follow, the top panel shows the data with the best-fit model superimposed. 
The lower panel shows the resulting x for each of the spectral channels, except in the high- 
resolution Chandra grating data where we plot the fiux residuals. 
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Fig. 2. — Same as in Figure 1 for time segment lb. 
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Fig. 3. — Same as in Figure 1 for GRB970828 time segments A and C 
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Fig. 4. — Same as in Figure 1 for GRB970508 time segment B. 
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Same as in Figure 1 for the dispersed spectrum of GRB991216. 
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Same as in Figure 1 for the zeroth order spectrum of GRB991216. 



> 



CO 
CI 

o 
o 



CO 

q 
d 

^^ 
o 



CM 
O 



O 
C\J 



GRB000210 
C/iandra/ACIS-S 



+\l. 



^^^^W 



H 1 1 h 



^++iHi+-Vt+;^+*+++^ 



02 

I 



0.5 1 2 

channel energy (keV) 

Fig. 7.— Same as in Figure 1 for GRB000210. 
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Fig. 8.— Same as in Figure 1 for GRB000214. 
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Fig. 9.— Same as in Figure 1 for GRB000926. 
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Fig. 10.— Same as in Figure 1 for GRB001025. 
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Fig. 11.— Same as in Figure 1 for GRB010220. 
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12.— Same as in Figure 1 for XRF011030. 
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Fig. 13. — Same as in Figure 1 for GRB011211 time segment 1. 
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Fig. 14. — Same as in Figure 1 for GRB011211 time segment 2. 
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Fig. 15.— Same as in Figure 1 for GRB020321. 
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16.— Same as in Figure 1 for GRB020322. 
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Fig. 17. — Same as in Figure 1 for the dispersed spectrum of GRB020405. 
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18. — Same as in Figure 1 for the zeroth order spectrum of GRB020405. 
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Fig. 19. — Same as in Figure 1 for the dispersed spectrum of GRB020813. 
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Same as in Figure 1 for the zeroth order spectrum of GRB020813 
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Fig. 22. — Same as in Figure 1 for the zeroth order spectrum of GRB021004. 
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Fig. 23.— Same as in Figure 1 for GRB030226. 
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24. — Same as in Figure 1 for the dispersed spectrum of GRB030227. 
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Fig. 25. — Same as in Figure 1 for the zeroth order spectrum of GRB030227. 
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Fig. 26. — Same as in Figure 1 for the dispersed spectrum of GRB030328. 
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Fig. 29.— Same as in Figure 1 for GRB030329 epoch 2. 
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Fig. 30. — Same as in Figure 1 for GRB031203 epoch 1. 
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Fig. 31. — Same as in Figure 1 for GRB031203 epoch 2. 
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Fig. 32.— Same as in Figure 1 for GRB040106. 
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34. — Same as in Figure 33 for GRB970508 section lb. 
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Fig. 35. — Same as in Figure 33 for GRB970828 GIS time section A and C 
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Fig. 36. — Same as in Figure 33 for GRB970828 SIS time section A and C 
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Fig. 37. — Same as in Figure 33 for GRB970828 GIS time section B. 
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Fig. 38.— Same as in Figure 33 for GRB970828 SIS time section B. 
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40.— Same as in Figure 33 for GRB991216 first order HEG. 
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41.— Same as in Figure 33 for GRB991216 first order MEG. 
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Fig. 42.— Same as in Figure 33 for GRB000210 ACIS. 
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Fig. 43.— Same as in Figure 33 for GRB000214 MEGS. 
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Fig. 44.— Same as in Figure 33 for GRB000926 ACIS. 
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Fig. 45.— Same as in Figure 33 for GRB001025 PN. 
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Fig. 46.— Same as in Figure 33 for GRB001025 MOS. 
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Fig. 47.— Same as in Figure 33 for GRBO 10220 PN. 
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Fig. 48.— Same as in Figure 33 for GRB011030 ACIS. 
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49.— Same as in Figure 33 for GRB011211 PN first 5 ksec. 
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Fig. 50.— Same as in Figure 33 for GRB011211 MOS first 5 ksec. 
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Fig. 51.— Same as in Figure 33 for GRB011211 PN t > 5 ksec. 
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Fig. 52.— Same as in Figure 33 for GRB011211 MOS t > 5 ksec. 
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Fig. 53.— Same as in Figure 33 for GRB020321 PN. 
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Fig. 54.— Same as in Figure 33 for GRB020321 MOS. 



> 


f\7 


v 




M 




m 




-P 








p 




o 




o 




O 


T-H 


T-t 




^'—^ 




^-^ 




U 




^Sm>^ 




O 





T 1 1 1 1 1 1 1 — [- 



— I 1 1 1 1 1 — I — r 

GRB020322 
XMM-Newton/Pli 




0.2 0.5 1 2 5 

channel energy (keV) 

Fig. 55.— Same as in Figure 33 for GRB020322 PN. 
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Fig. 56.— Same as in Figure 33 for GRB020322 MOS. 
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Fig. 57. — Same as in Figure 33 for GRB020405 zero order. 
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Fig. 58. — Same as in Figure 33 for GRB020405 dispersed. 
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Fig. 59. — Same as in Figure 33 for GRB020813 zero order. 
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Fig. 60. — Same as in Figure 33 for GRB020813 HEG dispersed spectrum. 
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Fig. 61. — Same as in Figure 33 for GRB020813 MEG dispersed spectrum. 
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Fig. 62. — Same as in Figure 33 for GRB021004 zero order. 
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Fig. 63. — Same as in Figure 33 for GRB021004 HEG dispersed spectrum. 
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Fig. 64. — Same as in Figure 33 for GRB021004 MEG dispersed spectrum. 
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Fig. 65.— Same as in Figure 33 for GRB030226 ACTS. 
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Fig. 66. — Same as in Figure 33 for GRB030227 PN cumulative spectrum. 
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Fig. 67. — A blow-up of Figure 66 in the 1.2 - 2.4 keV region. 
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Fig. 68. — A blow-up of Figure 66 in the 3.2 - 5.2 keV region. 
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Fig. 69. — Same as in Figure 33 for GRB030227 MOS cumulative spectrum. 
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Fig. 70. — A blow-up of Figure 69 in the 2.1-3.1 keV region. 
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Fig. 71.— Same as in Figure 33 for GRB030227 PN last 10.9 ksec. 
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Fig. 72.— Same as in Figure 33 for GRB30227 MOS last 10.9 ksec. 
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Fig. 73. — Same as in Figure 33 for GRB030328 zero order. 
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74.— Same as in Figure 33 for GRB030328 LEG dispersed. 
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Fig. 75. — Same as in Figure 33 for GRB030329 PN time segment 1. 
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Fig. 76. — Same as in Figure 33 for GRB030329 MOS time segment 1. 
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Fig. 77. — Same as in Figure 33 for GRB030329 PN time segment 2. 



-87- 





CO 




O 


^--v 




^H 




> 




(U 




^ 




m 


(M 


^ 


O 


S 




O 




o 




n 




O 




T-H 


^-1 


^-' 


o 


^n— S^ 




H 








O 





T 1 1 1 1 — r 



n 1 1 1 1 1 — r 



GRB030329 

XMM-Newton/UOS 

(epoch 2) 



/■ 



99.9% 

•v. 99% 

^ ^'\ data 




1 2 5 

channel energy (keV) 

Fig. 78.— Same as in Figure 33 for GRB030329 MOS time segment 2. 
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Fig. 79. — Same as in Figure 33 for GRB031203 PN time segment 1. 
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Fig. 80. — Same as in Figure 33 for GRB031203 MOS time segment 1. 
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81. — Same as in Figure 33 for GRB031203 PN time segment 2. 
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Fig. 82.— Same as in Figure 33 for GRB031203 MOS time segment 2. 
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83.— Same as in Figure 33 for GRB040106 PN. 
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Fig. 84.— Same as in Figure 33 for GRB040106 MOS. 



